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Abstract 

The stability theory of periodic traveling waves is much less advanced than for solitary waves, 
which were first studied by Boussinesq and have received a lot of attention in the last decades. In 
particular, despite recent breakthroughs regarding periodic waves in reaction-diffusion equations and 
viscous systems of conservation laws [Johnson-Noble-Rodrigues-Zumbrun, Invent math (2014)], the 
stability of periodic traveling wave solutions to dispersive PDEs with respect to ‘arbitrary’ pertur¬ 
bations is still widely open in the absence of a dissipation mechanism. The focus is put here on 
co-periodic stability of periodic waves, that is, stability with respect to perturbations of the same 
period as the wave, for KdV-like systems of one-dimensional Hamiltonian PDEs. Fairly general non- 
linearities are allowed in these systems, so as to include various models of mathematical physics, 
and this precludes complete integrability techniques. Stability criteria are derived and investigated 
first in a general abstract framework, and then applied to three basic examples that are very closely 
related, and ubiquitous in mathematical physics, namely, a quasilinear version of the generalized 
Korteweg-de Vries equation (qKdV), and the Euler-Korteweg system in both Eulerian coordinates 
(EKE) and in mass Lagrangian coordinates (EKL). Those criteria consist of a necessary condition for 
spectral stability, and of a sufficient condition for orbital stability. Both are expressed in terms of a 
single function, the abbreviated action integral along the orbits of waves in the phase plane, which is 
the counterpart of the solitary waves moment of instability introduced by Boussinesq. However, the 
resulting criteria are more complicated for periodic waves because they have more degrees of freedom 
than solitary waves, so that the action is a function of A -f 2 variables for a system of N PDEs, 
while the moment of instability is a function of the wave speed only once the endstate of the solitary 
wave is fixed. Regarding solitary waves, the celebrated Grillakis-Shatah-Strauss stability criteria 
amount to looking for the sign of the second derivative of the moment of instability with respect to 
the wave speed. For periodic waves, stability criteria involve all the second order, partial derivatives 
of the action. This had already been pointed out by various authors for some specific equations, 
in particular the generalized Korteweg-de Vries equation — which is special case of (qKdV) — but 
not from a general point of view, up to the authors’ knowledge. The most striking results obtained 
here can be summarized as: an odd value for the difference between N and the negative signature 
of the Hessian of the action implies spectral instability, whereas a negative signature of the same 
Hessian being equal to N implies orbital stability. Furthermore, it is shown that, when applied to the 
Euler-Korteweg system, this approach yields several interesting connexions between (EKE), (EKL), 
and (qKdV). More precisely, (EKE) and (EKL) share the same abbreviated action integral, which 
is related to that of (qKdV) in a simple way. This basically proves simultaneous stability in both 
formulations (EKE) and (EKL) — as one may reasonably expect from the physical point view —, 
which is interesting to know when these models are used for different phenomena — e.g. shallow 
water waves or nonlinear optics. In addition, stability in (EKE) and (EKL) is found to be linked 
to stability in the scalar equation (qKdV). Since the relevant stability criteria are merely encoded 
by the negative signature of (A -|- 2) x (A -|- 2) matrices, they can at least be checked numerically. 
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In practice, when Ai = 1 or 2, this can be done without even requiring an ODE solver. Various 
numerical experiments are presented, which clearly discriminate between stable cases and unstable 
cases for (qKdV), (EKE) and (EKL), thus confirming some known results for the generalized KdV 
equation and the Nonlinear Schrodinger equation, and pointing out some new results for more general 
(systems of) PDEs. 

Keywords: traveling wave, spectral stability, orbital stability, Hamiltonian dynamics, action, mass Lagrangian 
coordinates. 

AMS Subject Classifications: 35B10; 35B35; 35Q35; 35Q51; 35Q53; 37K05; 37K45. 
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1 Introduction 

Hamiltonian PDEs include a number of model equations in mathematical physics, like the (generalized) 
Korteweg-de Vries equation (KdV) or the Non-Linear Schrodinger equation (NLS). These equations and 
many others are known to admit rich families of planar traveling wave solutions, with more or less degrees 
of freedom. The most ‘rigid’ traveling waves are the so-called kinks, corresponding to heteroclinic orbits 
of the ODEs governing their profiles. Periodic traveling waves, which are the purpose of this paper, have 
the highest number of degrees of freedom. In between kinks and periodic waves in terms of degrees of 
freedom, we can find solitary waves, corresponding to homoclinic orbits. 

The actual existence of such waves follows from the Hamiltonian structure of the governing ODEs. We 
are most interested in their nonlinear stability, even though we can only hope for orbital stability, because 
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of translation invariance. The most efficient approach to tackle the orbital stability of Hamiltonian 
traveling waves has been known as the Grillakis-Shatah-Strauss (GSS) theory [GSS87], which provides a 
way of using a constrained energy as a Lyapunov function. This method crucially relies on the conservation 
of a quantity associated with translation invariance, termed impulse' by Benjamin [Ben84], and known 
as the momentum in the NLS literature. For solitary waves, the GSS theory provides a sufficient stability 
condition in terms of the convexity of the constrained energy as a function of the wave velocity. This 
constrained energy happens to correspond to what was called ^moment of instability^ by Boussinesq 
[Bou72] more than 140 years ago. Resurrected by Benjamin [Ben72] in the early ’70s, the ideas of 
Boussinesq have been made rigorous for many types of solitary waves in [GSS87, BSS87, BS88, BGDDJ05] 
(see also [BG13], [DBGRN14] and references in [AP09]). Together with the Evans functions techniques 
brought in by Pego and Weinstein [PW92], Kapitula and Sandstede [KS98], and many others, those pieces 
of work have led to a clear picture of which solitary waves are stable and which are not. 

By contrast, the theory is much less advanced regarding periodic waves. Apart from the higher number 
of degrees of freedom, the main difficulty comes from the fact that the nice variational framework set 
up by Grillakis, Shatah, and Strauss does not work for all kinds of perturbations of those waves. As a 
matter of fact, the theory of linear stability of periodic waves under ‘localized’ perturbations — that is, 
perturbations going to zero at infinity — is still in its infancy (see for instance [BD09, GH07b, BDNll] as 
regards spectral stability for KdV and the cubic NLS, and [Rodl5] for asymptotic linear stability of KdV 
waves), and the nonlinear stability under such perturbations is an open problem. In [BGNR14, BGNR13], 
the authors have contributed to the field by exhibiting several necessary conditions for the spectral 
stability of periodic waves in Hamiltonian PDEs. In particular, they have proved in a rather general 
setting that the hyperbolicity of the modulated equations ‘a la Whitham’ is necessary for the spectral 
stability of the underlying wave. More precisely, the existence of a nonreal eigenvalue for the modulated 
equations implies a sideband instability, which means that there are unstable modes for arbritrary small 
nonzero Floquet exponents. We shall not enter into details about these results here, we refer the reader 
to [BGNR13] and references therein — see also the recent related analysis in [JMMP14]. 

We are going to concentrate on the somehow easier problem of stability with respect to co-periodic 
perturbations, that is, perturbations of the same period as the wave (or equivalently, corresponding to a 
zero Floquet exponent). Our aim is to give as a clear picture as in the case of solitary waves, which by 
the way may be viewed as a limiting case of periodic waves — by letting their wavelength go to infinity. 
The great advantage of co-periodic perturbations is that they allow us to use the GSS approach in the 
simplest manner — by using basically only one additional conservation law (or constraint) to rule out 
the ‘bad’ directions from the variational framework —, and thus achieve nonlinear stability results at 
once^. This has been done for the cubic NLS by Gallay and Haragus [GH07a] — see also [BDNll, GP15] 
for more recent results, dealing with subharmonic perturbations, of which the period is a multiple of the 
period of the wave —, and for the generalized KdV (gKdV) by Johnson [Joh09] — for the classical or 
the modified KdV, see also [APBS06], [AP07, Arr09] for co-periodic, orbital stability and [DKIO], [DNll] 
for a more general result, which even handles subharmonic perturbations. Both NLS and gKdV can be 
viewed as specific cases of the abstract setting we are going to consider. Furthermore, this setting is built 
to include the Euler-Korteweg system, a fairly general model that is involved in various applications 
(superfluids, water waves, incompressible fluid dynamics, and nonlinear optics). The abstract systems we 
consider hold in one space dimension, and read 

9tU = /(Ejr[u]), (1) 

where the unknown U takes values in ^ is a skew-adjoint differential operator, and denotes the 
variational derivative of ^ Uj,) — the letter E standing for the Euler operator. In practice, we 

are most concerned with the case N = 2, which is the case for the various forms of the Euler-Korteweg 

^Even in the simplest context of co-periodic stability, it is indeed still much unclear how spectral stability can imply 
nonlinear stability in Hamiltonian frameworks. Concerning localized perturbations we even lack a clear notion of dispersive 
spectral stability that would be the analogue of diffusive spectral stability [Sch96, Sch98, JZIO, JNRZ14, RodlS]. 
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system, as well as NLS. In fact, (1) includes both the original formulation of NLS, with 



1 

0 


merely being the skew-symmetric matrix of Hamiltonian equations in ‘canonical’ coordinates, and its 
fluid formulation via the Madelung transform. In the latter case, 


^ = 'Rdx with B 



From now on, we assume that ^ = Jidx with B a symmetric and nonsingular matrix. This allows the 
case N = 1 with ^ = dx, which includes gKdV, and will enable us to make the connection with earlier 
results by Bronski, Johnson, and Kapitula. Furthermore, if iV = 2, we assume that the Hamiltonian Jtf 
splits as 

^ = jtf’iu] = jr(u, u^) = y{v, u) + c^{v, vx) , with u = 



and then that 


B'l 





Here above and throughout the paper, square brackets [•] signal a function of not only the dependent 
variable U but also of its derivatives \Jx, Vxx, •. • In this way, the abstract system in (1) reads as a 
system of conservation laws 

dtV = dx{BE.je[u]). ( 2 ) 

We recall that, when, as in cases under consideration in the present paper, Jif depends only on U and 
XJx, the Qf-th component (I < a < N) of the variational derivative EJ^[XJ] is 


(EJ^[U])„ 


dJif 

dUa 


(U,U,) - D, 


/ ajf 
\dU^,x 


(U,u, 


where D^, stands for the total derivative, so that 


D 


X 


/dJf_ 

\dUa,x 


(U,u, 


d^Ji^iV,Vx) 

dUisdUa,x 


Up^x 


+ 




Up 


,XX 5 


where we have used Einstein’s convention of summation over repeated indices. The main examples that 
fit the abstract framework in (2) are, besides the generalized Korteweg-de Vries equation 

(gKdV) dtv + dxp{v) = -dlv , 


and its quasilinear counterpart which is written in the more general form 

(qKdV) dtV = dx{E£[v ]), .e = £{v, Vx ), 
the Euler-Korteweg system in Eulerian coordinates, 

(EKE) I 

[ dtu-EudxU + dx{ES[p]) =0, = S{p,px), 

or in mass Lagrangian coordinates, 


(EKL) 


dsV = dyU, 


dsU = dy{E£:[v\) , £ = £{v,Vy). 
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We invite the reader to take a look at [BGNR13] for more details. 

The main results of the present paper are concerned with periodic traveling wave solutions to (2), 
with applications to (qKdV), (EKE) and (EKL). They consist of a sufficient condition for their orbital, 
co-periodic stability, and a necessary condition for their spectral, co-periodic stability. Both are expressed 
in terms of the Hessian of the constrained energy — to be defined in Section 2 hereafter — viewed as 
a function of the + 2 parameters determining periodic waves. The value of this constrained energy at 
a given wave profile happens to be interpreted as an abbreviated action integral along the corresponding 
orbit in the phase plane {(u,Ua;)}. Remarkably enough, as far as capillary fluids are concerned, that is 
for the systems (EKE) and (EKL) with an energy of the form 

S’ = F{p) + ]^jr{p)pl, £ = f{v) + ^K{y)vl, (3) 

the action integral physically corresponds to surface tension. The abbreviated action of periodic wave 
profiles also admits an interpretation in terms of the averaged equations for modulated wavetrains, in 
that it is dual to the wavenumber in the generalized Gibbs relation satisfied by the averaged energy (see 
Eqs (63) (64) in [BGNR14]) — this point of view is investigated further in a forthcoming paper. 

Co-periodic stability conditions are replacements for the — simpler — ones known for solitary waves. 
Indeed, while the abbreviated action depends on -I- 2 parameters for periodic waves, it depends on the 
sole solitary wave velocity once the endstate of solitary waves is fixed in in which case the abbreviated 
action merely coincides with the Boussinesq moment of instability. 

As is often the case, the necessary condition for spectral stability comes from an Evans function 
calculation, see Section 3. Phrased explicitely, it yields a sufficient condition for spectral instability, 
which is that the difference between N and the negative signature of the Hessian of the abbreviated 
action be odd. 

As to the sufficient condition for orbital stability, it relies on the GSS approach, together with a crucial 
algebraic relation, analogous to what has been pointed in [PSZ13] (see also [KP13, Proposition 5.3.1] and 
references therein). This relation makes the connection between the negative signatures of two sorts of 
Hessians associated with the constrained energy, namely, the differential operator obtained as the Hessian 
at the wave’s profile of the constrained energy viewed as a functional, and the {N -I-1) x {N + 1) matrix, 
corresponding to what Kapitula and Promislow [KP13] call a constraint matrix, and arising here as the 
Hessian of the abbreviated action integral under the constraint that the period of waves is fixed. This 
leads us to introduce a most important orbital stability index. All this is made more precise in Section 
2 below, and those necessary/sufficient stability criteria are actually derived in an abstract setting in 
Sections 3 and 4. 

Section 5 is then devoted to the application of these abstract results to (qKdV), (EKL), and (EKE), 
with energies as in (3). A striking result is that, in all these cases, our sufficient condition for orbital sta¬ 
bility mostly relies on the simple requirement that the negative signature of the Hessian of the abbreviated 
action be equal to N. 

In addition, we point out a close connexion between stability criteria for (qKdV) and for the Euler- 
Korteweg systems (EKL) and (EKE). We show in Section 5.2 that the systems (EKL) and (EKE) share 
the very same abbreviated action integral and constrained energy, in which the parameters of the waves 
turn out to be pairwise exchanged — as well as the constraints actually, if the period itself is considered 
as a constraint. This readily implies that our spectral stability criterion coincides for these systems. 
Furthermore, we prove that (EKE) and (EKL) actually have the same orbital stability index, equal to the 
negative signature of the Hessian of the abbreviated action minus two, even though the negative signatures 
of the unconstrained variational Hessians of respective Lagrangians can differ from each other. Regarding 
spectral stability with respect to ‘arbitrary’ perturbations — and not only co-periodic perturbations — 
we even show that the differential operators involved in the linearized systems associated with (EKE) and 
(EKL) are isospectral. Even though this seems a very natural result, the underlying conjugacy between 
eigenfunctions is far from being trivial. Moreover, we stress that the spectral conjugacy is not restricted 
to co-periodic boundary conditions and respects the Floquet exponent by Floquet exponent Bloch-wave 
decomposition. Both spectral and variational connections are pointed out here for the first time, up to 
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the authors’ knowledge. 

Some more specific examples — dictated by more or less classical choices of nonlinearities — are then 
investigated numerically in Section 5. This part relies very much on the fact that all our stability criteria 
are expressed in terms of the abbreviated action integral, which can be computed in the phase plane 
without any need of an ODE solver. Its derivatives are then computed by means of hnite differences. The 
coexistence of two grids of discretization — one for the integral and one for finite differences, and the high 
condition number of the Hessian matrices that are to be computed, induce some numerical difficulties 
that have been coped with by a suitable choice of mesh sizes. Numerous numerical experiments have 
been conducted, and their results are in accordance with those that can be computed analytically. In 
particular, our routine for computing the Hessian of the abbreviated action integral enables us to recover 
in a very precise way — and up to the small amplitude limit and to the soliton limit — the eigenvalues 
of modulated equations associated with some well-known completely integrable PDEs (namely, KdV, 
mKdV, and cubic NTS), as displayed in a forthcoming paper [BGMR15]. 

Coming back to the analytical part of our results, let us mention the following, important difficulty. 
In order to actually prove some orbital stability, we need a suitable local well-posedness theory for the 
Cauchy problem. This kind of result is of course heavily model-dependent. If there is for instance a huge 
literature on (g)KdV, it does not seem that anyone ever looked at the Cauchy problem for (qKdV) when 
the ‘capillarity’ factor k in 

f{v) + ^k{v)vI 

is not constant. This is done in a forthcoming paper [Miel5]. Regarding (EKE) and (EKL), still with 
energies as in (3) with variable and k, what we need is a basic adaptation to the ID torus of earlier 
results dealing with the Cauchy problem on the whole real line [BCDD06, BCDD07]. 

2 Summary of main results 

In order to define the constrained energy, we observe that the system (2) formally admits the additional 
conservation law 

a*^(U)=5,(J5"[U]) (4) 

with 

i?(U) := i U • B-^U, yiV] := U • + LJf[V] , 

L^[U] := . 

oc,x 

The dots • in the definitions of i? and y’ are for the ‘canonical’ inner product U • V = [/qEq, in 
Recall that the sans-serif letter E stands for the Euler operator defining variational derivatives. As to 
the other sans-serif letter L, it stands for a crude version — without any change of variables involved — 
of the "Legendre transform’, just defined by the formula above. We see that is associated with the 
invariance of (2) with respect to a;-translations because of the algebraic relation 

9,U = 9,(BEi?[U]). (5) 

As a consequence, for a travelling wave U = U(a; — ct) of speed c to be solution to (2), one must have by 
(5) that 

a,(E(jr + c^)[u]) = o, 

or equivalently, there must exist A G such that 

E(jr-Ec^)[U] = A. (6) 

Equation (6) is the Euler-Lagrange equation associated with the Lagrangian 

^ U,,; A, c) := U,,) -E c^(U) - A • U . 
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From place to place we shall refer to the components of A as Lagrange multipliers. We thus see that 
is a first integral of the profile ODEs in (6). We also notice that (6) implies 

L^[U] = ^[U] + c^(U), 

which is of course consistent with the conservation law in (4). In this way, the possible profiles U are 
determined by the equations in (6) together with 

Lif[U]=/x, (7) 

where p is a, constant of integration, which we shall sometimes refer to as an energy level, since is 
the conserved ‘energy’ associated with the Lagrangian L^. 

All this roughly shows that a periodic travelling profile U is ‘generically’ parametrized by {p, A, c) G 
]g7V+2_ Furthermore, if such a profile U is of period S, we can define the constrained energy of all 
S-periodic, smooth enough functions U by 

^[\5-,p,X,c] := [ (^^(U, U,,) + c^(U) - A • U + Ai) dx. 

Jo 

Denoting by 

0(/i, A, c) := ^[U;/r, A, c] = f u) + ci?(U) — A • U +/i) dx , (8) 

Jo 

we can see by using (7) and a straightforward change of variable that 

/ rte 

— (u,u,,)dv (9) 

is the abbreviated action of S along the orbit described by v in the (u, Ua;)-plane. This is the reason why, 
as was observed in [BGNR13, Proposition 1], the partial derivatives of 0 are merely given by 

0^ = S,Va 0 = -/ Udx,0c=/ i2(U)da;. (10) 

Jo Jo 

We can now state our stability conditions in a more precise way. 

Necessary condition for spectral, co-periodic stability. A periodic traveling wave solution to 
(1), of profile U and period S, is said spectrally stable with respect to co-periodic perturbations if the 
spectrum of the linearized operator 


.!^:= /Hess(.^-hc^)[U] 

in L^(IR/E!Z) is purely imaginary. We cannot expect more than this neutral stability, because of symme¬ 
tries^. A necessary condition for spectral, co-periodic stability is 

• det(Hess0) < 0 in the case A^ = 1, 

• det(Hess0) > 0 in the case N = 2. 

The scalar case A^ = 1 is a generalization to quasilinear equations of the results found by Bronski and 
Johnson [BJIO]. For both cases, the proof is based on the fact that possible unstable eigenvalues z are 
characterized by D{z) = 0, with 

D{z) :=det(F(S;z)-F(0;z)) 

^The eigenvalues of ^ outside real and imaginary axes arise as quadruplets (. 2 , 2 ,— 2 ,— 2 ), nonzero real or imaginary 
eigenvalues come in pairs (. 2 , —z). 
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where F(-; z) denotes the fundamental solution of the ODEs in £/U = zlJ, and on asymptotic expansions 
showing that when r is real 

D(r) (—I det(B~^)| det(Hess6) + and D{r) > 0 for r » 1. 

These results are collected in a rigorous manner in Theorem 1 (for N = 2) and Theorem 2 (for N = 1) 
in Section 3. 


Sufficient condition for orbital, co-periodic stability. It is obtained through a variational argu¬ 
ment. We assume that ^ 0, and dehne the constraint matrix as 


C := 


V0^(8) ve^ 

©MM 




with V being a shortcut for the gradient with respect to (A, c) at fixed /x. (Note that the coefficients of 
C are made, up to a factor —1/0^^, of all the 2x2 minors of V^0 containing 0^^.) Let us denote by 
A the differential operator obtained as the Hessian of the constrained Hamiltonian: 


A := Hess(Jf’-hci?)[U]. 

If C is nonsingular, and if the negative signatures of the operator A and of the matrix C happen to 
coincide, then the periodic travelling wave solutions to (2) of profile U are orbitally stable in T^(IR/SZ). 
This is the purpose of Theorem 3 in Section 4. Its proof is based on the following algebraic relation, 
shown in [BGNR13] (see also [PSZ13], [KP13, Proposition 5.3.1] for similar observations), 

n(A) = n(A|Tu<^) + n(C), 

where n denotes negative signature, and is the tangent subspace to the {N + 1) codimensional 

manifold 

{U e l2(M/SZ) ; /“ ^(U) dx = ^(U) da;, fg" U dx = /q" Udx} . 

(The space Tjj^ actually corresponds to what Kapitula and Promislow [KP13] call an admissible space.) 
According to that relation between negative signatures, the fact that n(A) = n(C) implies that the 
operator A\Xu‘^ i® nonnegative, and up to factoring out the null direction TJ_^, this roughly shows that 
the functional n, A, c] has a local minimum at U and its translates U(- -I- s) on 

{U e l 2(M/SZ) ; /“ ^(U) dx = ^(U) dx, U dx = /q" Udx} . 

Orbital stability can then be achieved by a contradiction argument as in [BSS87, GSS87], see [BGNR13], 
or in a direct way as in [Joh09, GH07b] (see also [DBGRN14, §4.2 & §7.3] for an interesting discussion of 
the pros and cons of these arguments). We choose the direct way for the proof of Theorem 3 in Section 

4. 

In practice, we need to evaluate the orbital stability index n(A) — n(C). For the ‘concrete’ systems we 
consider in Section 5, we can infer from a Sturm-Liouville argument that n(A) € {1,2}. In particular, 
extending results by Johnson [Joh09] for (gKV), we show that 0^^ > 0 implies n(A) = 1. Furthermore, 
we see that when N = 1, 0^^ det(C) = det(Hess0). Hence a more explicit — but partial — version of 
the sufficient stability condition: 

©Ai/i > 0, det(Hess0) < 0, 

which ensures indeed that n(A) = n(C), and is of course consistent with Johnson’s findings in the 
semilinear case. In fact, we recover in Section 5.1 the more general sufficient condition for n(A) = n(C) 
in (gKdV) that was later derived by Bronski, Johnson and Kapitula [BJKll] — using in particular that 
Qfiij. < 0 implies n(A) = 2 —, namely 


©MM 0 • 


©AA ©mA 
©A/i ©MM 


7 ^ 0, det(Hess0) ^ 0, n(Hess0) = 1. 



More generally, for our concrete systems, we prove by related arguments — Sturm-Liouville theory (see 
Lemma 2) and simple algebraic relations (see Proposition 1) — the remarkable identity for the orbital 
stability index 

n(A) — n(C) = n(Hess0) — N. 

Therefore, our nonlinear stability result applies as soon as 0^^ ^ 0, det(Hess0) ^ 0 and 

n(Hess0) = N. 

Our approach is applied to (qKdV) in Section 5.1, Theorem 5. 

Details regarding the systems which have motivated this work, (EKL) and (EKE), are given in Section 
5.2. A most important fact is that these systems share the very same abbreviated action, 0(/i, A,j, ct) 
defined in (35), where a is the speed of EKE waves, j is the ‘speed’^ of EKL waves, and the roles of the 
parameters ^ and A are exchanged when we go from (EKL) to (EKE) and vice versa: fj, is an energy level 
for EKE waves and a Lagrange multiplier for EKL waves, A is a Lagrange multiplier for EKE waves and 
an energy level for EKL waves^. As a consequence, the stability criteria which are expressed only in terms 
of det(Hess0) and n(Hess0) coincide for corresponding EKE waves and EKL waves. By contrast, the 
individual negative signatures n(A) and n(C) are in general not preserved by going from one formulation 
to the other. This explains why some simplified, partial criteria — analogous to those in [Joh09] for 
(gKdV) for instance — are actually formulation-dependent. 

The fact that waves should be simultaneously stable in both formulations seems very natural from 
a physical point of view. However, this is not that obvious to prove mathematically, because the mass 
Lagrangian coordinates are obtained from the Eulerian coordinates through a nonlinear and nonlocal 
change of variables that depends on the solution itself. 

As far as spectral stability is concerned, and not only co-periodic stability actually, we prove in 
Theorem 6 (also see Remark 6) that the corresponding linearized operators are indeed isospectral. The 
proof is quite simple once we reformulate the nonlinear problems in a suitable way, but it reveals that 
the kind of conjugacy between those operators is not trivial. We are not aware of any earlier result of 
this type. 

As regards the co-periodic orbital stability, its simultaneous occurrence in both formulations (EKL) 
and (EKE) is supported by the idea that corresponding EKE waves and EKL waves share the same 
constrained functional pt, X,j,a], and that the constraints are preserved by passing from Eulerian 
coordinates to mass Lagrangian coordinates^. If the vanishing of our orbital stability index were exactly 
equivalent to the fact that the functional pi, A, j, cr] has a local minimum on the constrained manifold 

at the wave profile and its translates, this would show the equivalence of co-periodic orbital stability 
for (EKL) and (EKE). We find out that the issue is a little more subtle by looking at our abstract result. 
Theorem 3. Recalling that the roles of the ‘concrete’ parameters pi and A are exchanged when we go from 
(EKE) to (EKL), we see that the main assumptions for applying Theorem 3 to (EKE) and (EKL) are, 

0Ai/i 7^ 0, det(Hess0) ^ 0, n(Hess0) = 2, 

for the former (see Theorem 11), and 

0AA 7 ^ 0, det(Hess0) ^ 0, n(Hess0) = 2 

for the latter (see Theorem 10). The slight discrepancy between these two sets of conditions obviously 
comes from the derivatives 0^^^ and 0 aa, which correspond respectively to the derivative with respect 
to /i of the wave period in Eulerian coordinates — \i being the energy level in these coordinates — and 

®We use some quotes here because this is not a speed from the physical point of view, it actually corresponds to a mass 
transfer flux across the corresponding EKE waves. 

■^To avoid the introduction of too many notations, we have chosen to use the greek letters p, and A with a meaning in the 
‘concrete’ examples (EKL) and (EKE) that is slightly different from their meaning in the abstract framework, see Table 2. 

®In fact, this is true provided that we also consider the period as a constraint, and thus prescribe the (N + 2) constraints 
fg dx = S, fg U(x)d3:: = /(f U(a;)da:, fg ^(U(a;))dx = f^^(Uix)dx. 
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the derivative with respect to A of the wave period in mass Lagrangian coordinates — A being the energy 
level in these coordinates. They are not a priori related to each other. However, as long as we regard the 
vanishing of either 0^^^ or Qx\ as anomalous transitions, we can indeed think of the periodic waves as 
being simultaneously orbitally stable with respect to co-periodic perturbations in Eulerian coordinates 
and mass Lagrangian coordinates. Again, this is not that an obvious result, because the meaning of 
co-periodic perturbations is different from one formulation to the other®. 

Finally, it turns out from a simple algebraic computation that for (EKL) the negative signature of 
the Hessian of the constrained Hamiltonian in mass Lagrangian coordinates coincides with the negative 
signature of the qKdV operator, 

a := Hess(£-|- cxj)[y\ , 

where xj := is the qKdV impulse, and c = —is prescribed by the speed of the EKL wave. Two 
ingredients then lead to a set of sufficient stability conditions for (EKL). The first one is that, as mentioned 
above, n(a) is known in terms of the sign of 0^^^, where 9 is an alternative notation for the abbreviated 
action associated with qKdV waves — to avoid confusion with the one associated with EK waves, still 
denoted by 0. The second ingredient follows from the observation that 0 is explicitly related to 0, so 
that by the chain rule Hess0 can be expressed in terms of Hess0 and V0. 


3 Co-periodic, spectral instability 

3.1 General setting 

As in the introduction, we consider an abstract Hamiltonian system as in (1), with 


U = 


V 

u 


/ = B5, , B-i 



5 7 ^ 0 , 


jr = jr[u] = ^^(u, u,,) = j{v, u) + ^{v, v ^), 
and denote by ^ the impulse — or momentum — defined by 

^(U) = iu-B-^U. 

Furthermore, let us assume that is strongly convex with respect to u, and that S’ is strongly convex 
with respect to Vx- Both (EKE) and (EKL) fit this abstract framework, with 


U = ^ ^ , J^(p,w) = , <S'{p,Px) = ^^{v)pI + F{p) 

for the former, and 

U = ^ ^ , y{v, w) = (u, Vy) = ]^k{v)vI + f{v) 


for the latter, as long as k, Jf, and p take positive values. 

Recall that profiles U of periodic wave solutions to (1) are characterized by the algebro-differential 
system made of (6), (7), which depends on the parameters (/x. A, c) € Equivalently, by using (5), we 
can view the prohle U of waves of speed c as a spatially periodic, and steady solution to the system (1) 
rewritten in a mobile frame, which reads 


dtV = B dx{E{J^ + c^)[U]). 


( 11 ) 


®More precisely, the prescription of the period on one side corresponds to a ‘zero mass’ perturbation on the other side, 
see Remark 9 for more details. 
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For later use, let us note that (11) admits the conservation law 


dt^iU) = d^{U ■ + c^)[U] + + c^)[\J]), (12) 

which is just (4) written in the mobile frame. In a similar way, let us note that as soon as Eq. (6) holds 
true, Eq. (7) equivalently reads 

U-E(Jf’ + ci?)[U] + l{J^ + c^)[U] = Ai, (13) 


We now fix such a periodic profile U, say of period S, and assume without loss of generality that 
vanishes at a; = 0, which will simplify a little bit our computations. Linearizing (11) about U, we receive 
the following system, in which the same notation U now stands for the variation of the original U around 

U, 

5tU = Ba,(AU), (14) 


where 

A := Hess(^ + c^)[U]. 

In general, HessJ^[U] is the differential operator defined by 


(Hess^[U]U)a = 


dUcdUp 


U0 


a2^[u] 


UjH.x 


a2^[u] 


Up + 


dUa,xdUp^, 


■Up,. 


and similarly for However, Hessi?[U] happens to merely coincide here with the matrix and A 

is like HessJ^[U] a self-adjoint differential operator, of second order in v, with periodic coefficients of 
period S. Our main purpose here is to derive a criterion ensuring that the composite differential operator 


:= B i9a;A 


does not have any spectrum in the complex, right-half plane, when the eigenfunctions are sought El- 
periodic. 

Let us recall the classical observations that the profile equation in (6) implies, by differentiation in x, 
that AU^ = 0, and by differentiation with respect to parameters /r, A, c, we see that^ 

= 0, =0, .c/U;,, = 0, .2/U, = -U, . 

Another useful relation, which follows from (14) but is more easily derived by linearizing (12) about 
U, is the following 

dt(V^(U) ■ U) = d^(U ■ E(,J^ -f ci?)[U] -E U • AU - U • Vu(^ + c^)[U]) 

+ da; Va.V + dl^S[v\ Vo- V^;) , 

which can be simplified into 

(9t(V^(U) • U) = a^(U • AU - KVa;xV + AH,, u^) (15) 

where k := d‘l^(^[y\ > 0 by assumption. 

For all z € C, let us consider the spectral problem associated with (14), 

zU = (16) 

which amounts to looking for solutions of (14) of the form e^‘U(a;). For such a solution, (15) readily 
implies by integration the following relation 


z f Vi2(U(a:)) • U(a:) dx = U(0) • [AU] - ^(0) u,,(0) [u], 
^0 


^Throughout the paper, subscripts fi, Ai, A 2 — or simply A when = 1 —, and c stand for partial derivatives with 
respect to the parameters p, A, c. 
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where we have used the shortcut [/] for expressions of the form /(S) — /(O), and the fact that U is 
S-periodic and chosen such that ^^.(O) = 0. Furthermore, by also integrating (16) over a period, we 
obtain 

z [ U(x) dx = B [AU]. (17) 

Jo 


Therefore, we have 


z [ V^(U(x)) • U(x) dx = z B-^iU(O) • [ U(x) dx - k(0) v^^{0) [x] . 
^0 Jo 


This relation will be used in a crucial way below, given that ^(0) > 0 and ^ 0. 


(18) 


3.2 Evans function 

The eigenvalue equations in (16) consist of a system of ODEs, which is of first order in u and third 
order in v. Rewriting this system as a first-order system of four ODEs — linear ODEs with S-periodic 
coefficients, and denoting by F(-;z) its fundamental solution, we see that the existence of a S-periodic, 
nontrivial solution to (16) is equivalent to D{z) = 0, where 

D(z) := det(F(S;z)-F(0;z)). 

The function D : C —>■ C is called an Evans function®. 

Theorem 1. In the framework of Sections 2 and 3.1, we make the ‘generic’ assumption 

(HO) There exists an open set D and a family of periodic traveling profiles U smoothly parametrized 

by {pL,\,c) G such that (6)-(7) hold true, that is 

E(Jf' + c^)[U] = X, U • Ejr[U] -b L.^[U] -b ci?(U) = /x. 


If N = 2 then the Evans function D defined above has the following asymptotic behaviors 

D{r) r^ (— det(B”^)) det(Hess0) -b o(r^), D{r) > 0 for r 1. (19) 

//det(Hess0) < 0, then the corresponding wave is spectrally unstable. 

Notice that det(B^^) = —so that (—det(B~^)) = |det(B^^)|. 

Proof. We begin by observing that U^,, U^, U;),^, are solutions to the ODEs in (16) with z = 0, as a 
consequence — by differentiation in x — of the relations 

AU, = 0, AU^ = 0, AU;,^ = (?))’ = ( 1 ) ’ 

which themselves come from the differentiation of (6). Furthermore, (U,, U^, is an indepen¬ 
dent family. Indeed, would a linear combination /3oU, -b /daU^ + fiiHx., + zero, the last two 

equations in (20) above would imply fii = = 0, so that /^oU, -b /^aU^ = 0, which turns out to be 

impossible unless fio = fis = 0- Indeed, by differentiation of (7) we see that 


-fSVxxlfu + -EVxxTx + EVxV.xx = 0> 


so that and x, cannot be colinear. (In the computation here above, we actually have differentiated 
(7) under the form (13), and used the same simplification as in the derivation of (15), as well as the first 

®Here with Floquet exponent equal to zero, since we search for co-periodic eigenfunctions only; a more general Evans 
function would be D{z, a) := det(F(S; z) — e'“F(0; z)), defined for all Floquet exponents a £ M. 
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two equations in (20) to cancel out the terms U • AU^, and IJ • AU^, even though these simplifications 
are not necessary to show that and cannot be colinear. If the second equation obtained in this 
way is trivial, this is not the case for the first one.) Since we have enforced ^^.(O) = 0, the above relations 
imply in particular 

f£(0)u,,,,(0)u^(0) = -1. (21) 

This preliminary observation that (U,j,,U^,U_,^^,U;^^) is an independent family allows us to consider 
the family of independent solutions z), j = 1, 2, 3,4, to (16) defined by the initial conditions 


{v^ (0; z), vi{0\ z), z), (0; z))"^ 


(Vx(0),Vxx(0),llxa:^(0),M,,(0))'r if j = 1, 

if j = 2, 

(llAi(0)>llAi.x(0)>l'Ai,a:a:(0)>MAi(0))'^ if J=3, 

(llA2(0)>llA2.x(0)>l'A2.a:^(0)>llA2(0))^ if i=4- 


We see that the Evans function equivalently reads D(z) = E{z)/A, where 

n^O) ••• u4(0) 

v^(0) ••• u;^(0) 

^'L(o) ••• i’L(o) 

1(1(0) ••• u4(0) 

Here above and in the remaining part of this proof, the notation [•] is reserved for differences of values 
between a; = S and a; = 0, and has nothing to do with the evaluation of functionals anymore. 


E{z) := 


.,,11 

r .,,41 

V \ 

• [ V \ 

ll 

r 4i 

•• 

• [ Kl 

.,,1 1 

r 4 1 


• ) 

ll 

r 4i 

It J 

• [ 


A 


Low frequency expansion. This is basically a variation on the computation made in [BGNR14, 
Section B.2], with a few more details for the reader’s convenience. We begin by observing that 


[AU] 


-k(0) [Vxx] \ * * * 

0(0) [m] y V * 0 0 



where a := d^J^{v,u) > 0 by assumption, and * stand for immaterial real numbers coming from the eval¬ 
uation of U and its derivatives at 0. Therefore, we can make some row combinations in the determinant 
defining E{z) by using (17), and obtain 


E{z) = —z^det(B 1 )k( 0) ^ a(0) ^ 




We can proceed in a similar way by using (18), which yields 


E{z) = z3det(B-i)«(0)-2a(0)-iu,,(0)-i 

/,^V^(U)-Ui ... /“V^(U).U4 

X ••• 

/o"ui ... /“U^ 

Now, we observe that (Ui)p^o = Us? and by differentiating zTJi = we see that U,,, = 

£/(dz'U^)lz^o. Therefore, IJi -f is in the kernel of jzf, which is spanned by (U^,, U^, U;,^^, U;^^) 
since this is an independent family of solutions to ^/U = 0, which is equivalent to a first-order system of 
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four ODEs. We also have that (U^)| 2 =^o = U;i> (U^)| 2 =o = Uai) (U^)| 2 =o = by the choice of initial 
conditions. We thus find that 

E{z) = -z4det(B-i)«(0)-2a(0)-i^;,,(0)-i 

/“V^(U).U, /,^V^(U).U^ /“V^(U).U,^ /„-Vi2(U).U,^ 

X [V-cJ h^i,x] [tlAi.x] hx^J +0{z^), 

rn, /“U,^ 


which equivalently reads, according to (10), 


E(z) = —det(B ^)k(O) ^a(O) ^ ^ x 


©cc 

©cp 

BcAi 

BcA2 

[—c,ai] 

[—/j.,ai] 

fcAi,rr] 

bA2,x] 

BAic 

Bai/j. 

Ai Ai 

Bai A 2 

Ba2C 

Ba2/J- 

Ba 2 Ai 

Ba2 A 2 


+ 0{z^). 


Finally, using that [v^ = —^a'Hxxi^)^ for a = c, /r, Ai, A 2 — which merely comes from the differentiation 

of the relation v^{^) = tta;(0) with respect to those parameters — we obtain 

^CC ^Cfl ©cAi ®cA2 

E(z) =z4det(B-i)K(0)-2a(0)-i X ^ Xi a^ 

Oaic f^AiAi OaiA2 

f^A2C Oa 2^ f^A2Ai OA 2 A 2 

that is, recalling also from (10) that S = 0^, 

E{z) = z"* det(B^^) k(0)“^ q;( 0)”^ det Hess© + 0{z^). 


On the other hand, using that 
AU(0) = 


-k(0)u2,2,(0) \ * * * 

a(o)u(o) y V * 0 0 


u(0) \ 

Wa;(0) 

uiO) ) 


in general, and Eqs in (20), we can compute explicitly A by making row combinations again. By using 
(21) in the last step, we thus find that 


lla;(0) 


HaAO) 

lfA2(0) 

VxxiO) 

Vt^,xi0) 

lfAi,a;(0) 

11A2.x(0) 

—XXX (f^) 

'!lf,,xxi^) 

llAi.xa:(0) 

i'X2,xx(0) 

Ma;(0) 

u^(0) 

Mai(0) 

Ma2(0) 


0 11^(0) tlA,(0) lfA2(0) 

_ —a:a:(f^) —/^,a:(f^) —Ai,a:(f^) —A2,a:(f^) 

tiiuj ayy)} 0 0 10 

0 0 0 1 
= K(0)"^a(0)-i u,,,^(0)^^(0) 


= “15(0) ^a(O) 
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Altogether with the expansion of E{z), this gives 


D{z) -04det(B-i) det Hesse + 0{z^) . 


High frequency expansion. In order to find the sign of D{z) for z G M+ large enough, we can invoke 
a homotopy argument^ . For this argument to work out, we must check that there exists i? > 0 so that 
D{r) 7 ^ 0 for r > R, and moreover that this R can be found to be uniform along the family of Evans 
functions associated with a continuous path going from the operator to a, simpler operator, say saf, for 
which we can compute the Evans function D explicitly. We choose 


£/ :=Bd^ Hess.^[U], 



^[U] := 1^2 + 



and postpone the search for R to the next paragraph. 

Once we have this i?, we know that D{r) and D{r) have the same sign for r > R. Let us compute 
D{r). The eigenvalue equations jz/U = z\J equivalently read 

J —d^v = zbu, 

\ dxU = zbv, 


or 


/ V \ 

Vx 

'^xx 

\ u J 


( y \ 

Vx 

Vxx 

\ U J 


A{z) 


/ 0 1 0 0 \ 

0 0 1 0 

0 0 0 

\ zb 0 0 0 J 


It is then a simple exercise to compute D{z) = det(e“^(^^ ~ !)• In particular, for z = 
has four distinct eigenvalues, ±(1 ± i)y^^r\b\j2. Therefore, it is diagonalizable and 


r e (0, +oo), A(r) 


det(e“'^^’'^ — 1) = (e“ ’’l'>l/2 _ l) (gS {l-i)y/r\b\/2 _ 2 ^^ {-l+i)y/r\b\/2 _ ^gB (-l-i)y/r\b\/2 _ 


gB(l+i)v/r|&|/2 


1 gB(-l+i)v/r|&|/2 


2 

1 > 0 . 


High enough frequencies are not eigenvalues. For any 9 € [0, 1], we set 

:= Be dx Hess.^[U] where JFe = OM’+ {1 - 9 )jF and Bg^ = 6iB-i + (1 - e)B-\ 


Notice that this does define a B^ since the formula for B^^ defines a matrix with determinant det(B~^). 
The aim is to find R G (0, +oo) such that, for all 9 G [0,1], the operator £/e does not have any 5-periodic 
eigenfunction associated with a real eigenvalue r > R. Observing that is — on purpose — exactly of 
the same form as .e/, by just replacing B’ -1, .y, ^ by 



we can drop 9 and seek R for under the only assumptions that b is fixed and positive, a may vary 
while staying bounded, k = and a = dinv , u] are bounded, positive and bounded away from 

zero. Such an R can be derived from some rough a priori estimates. A similar computation was made in 
[BGNR14, Section B.l], with a slight mistake which can be fixed by modifying the high order estimate 
accordingly with what follows. 

®Usual in many related computations, but apparently used for the first time in a periodic context ; see Remark 2. 
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Let us write the eigenvalue equations (16) in a more explicit way, 

f i9a;(Mn + + yrt) + cdx{av + bu) = r{av + bu), 

\ dx{':iv + au) + c dx{bv) = r bv , 

where M := Hess<f[:c] = —dxKdx + S with g bounded, and ^ := dxJ^[v,u\, j := duv-^[v,y] bounded too. 
On the one hand, taking the inner product of the system above with 

b(") = ( , ,0 

\ u J y v/b— au/b J 

in L^(M/SZ) and integrating by parts, we find that 

^(Ikllis + IklliO = Jo ( Uxdxinvx) + udx((g +/3)v + :fu) + (v - au) dxijv + au))/b dx 

for some constant C depending only the bounds on a, k, g, a, /3, j. On the other hand, taking the inner 
product of the system above with (vx,Ux)'’^, and integrating by parts again, we obtain 

0 = /g“ ( Vxx dx(KVx) + Vx dx((g + /3)v + + Ux dxi^v + au) 

+cvx dx(av + bu) + cux dx(bv)) dx 

= Jq {kvIx + aul) dx + /g“ ( KxVxVxx -IVxxU +Vx dxiig + 13 + ac)v) 

+ Uxdx{{'y + 26c) x) dx 

^ I fo i'^^xx + aul) da; - C (||x||^i + \\u\\\2) 

for some other constant C depending only the bounds on a, k, g, a, 13, 7 , 7 ^. Using once more that 
K, and a are positive and bounded away from zero, we thus find a constant C" such that 

\\vrH. + \\urm < c"{\\v\\h + \\u\\h). 

(Note that ||x 2,||^2 has been absorbed in the left-hand side.) Therefore, we have 

KIMIi= + lhlli.)<c^c'"(||x||i 2 + ikiii.), 

which implies that v and u must be zero if r > CC". 

Conclusion. By the mean value theorem, a necessary condition for stability of the wave is that D does 
not change sign on (0, -foo). Combining the low frequency expansion with the fact that Dir) is positive 
for large r, we obtain the necessary condition for co-periodic stability 

— det(B~^) det Hess© > 0. 

which requires that det Hess© be nonnegative. □ 

We can show a similar result in the case iV = 1, which corresponds to (qKdV), or equivalently 
U = X, ,y^\y+\ = e.{v,Vx), B = 1, ^ = £j = in the abstract form (2). In this case, the profile 
equation ( 6 ) reduces to the scalar Euler-Lagrange equation 

E{£+cxj)[v\ = A, 

and the abbreviated action is defined as in ( 8 ) by just dropping the u-component and 

fT 

9{^,X,c):= J {£{v,Vx) + cx^v) — Xv +fi) dx . 
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For a reason that will be clarified in Section 5.2, we have substituted T for as the period of the wave. 
The Evans function is also defined as before by 


d{z) := det(f(T; z) — f(0; 2 )), 

where is now the fundamental solution of the third order ODE 

(9a;(Hess(.e+ cxj)[y\v) = z v 


( 22 ) 


viewed as a first-order system. 

Theorem 2. Under assumption (HO), in the case N = 1, the Evans function defined above has the 
following asymptotic behaviors 

d{r) r^ det(Hess0) -|- o(r^), d{r) < 0 for r » 1. (23) 

//det(Hess0) > 0, then the corresponding wave is spectrally unstable. 

For consistency, notice that here det(B“^) = 1. 

Proof. We can basically copy-paste computations from the proof of Theorem 1, by taking B = 1, and 
dropping the u-components, J^, and all terms in A 2 . We just have to pay attention to where the signs 
change. A first, obvious one is that det B is now positive. There is another change of sign in the 
computation of E{z) near zero, because there is now only one row where we find a minus sign by writing 
f^vdx = —0\. So these two changes of sign give the claimed asymptotic expansion at zero. There is a 
third, and last change of sign in the computation of d{r) for large r. Indeed, the ODE to solve is now 

—d^v = rv, 

which has the three wavenumbers fco = ki = and ki for r > 0. Therefore, d{r) = 

(e'=« - 1) le'^i - Ip <0. □ 

Remark 1. Our main implicit restriction — even at the abstract level — is that we only consider 
systems for which, by a suitable number of integrations, the original traveling-wave profile system may 
be converted in a planar Hamiltonian, reduced, profile equation. Otherwise one would expect, as in the 
well-studied — and algebraically much simpler — case of solitary waves, some orientation index to enter 
in formulas for stability indices. For solitary waves, the computation of the necessary orientation index 
from geometric invariants is still the object of intense research ; see for instance [CB15] and references 
therein. 

Remark 2. As said in the introduction, this result confirms earlier findings by Johnson [Joh09] in 
the special case when k is constant. Though our presentation for (qKdV) does not follow the one 
by Johnson for (gKdV), some of our steps do not differ significantly. By contrast, some others are 
fundamentally different, as required by the quasilinear nature of our problem. For instance, in the 
foregoing proof, the purpose of the homotopy argument and the auxiliary resolvent estimates is precisely 
to reduce computations to a semilinear case. When equations are already in semilinear form, those 
techniques are not needed, and a readily regular limit A —>■ 00 leads to a constant-coefficient problem. 
Likewise, the local well-posedness results invoked in Section 5 for applying our abstract orbital stability 
to actual PDFs is dramatically improved for semilinear versions of those equations. These observations 
are instances of the usual rule of thumb that departures of quasilinear strategies from semilinear ones are 
only required when some high-frequency control is needed. 

Remark 3. It is instructive to seek parallels of our results in the classical stability theory for steady 
states of finite-dimensional Hamiltonian systems of ordinary differential equations. Indeed, up to replacing 
Evans’ functions with characteristic polynomials, the foregoing proofs echo the classical proof that steady 
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states at which the Hessian of the Hamiltonian is nonsingular and has an odd number of negative directions 
are spectrally instable. We claim that the analogy goes further. On the one hand, it follows from our proof 
that the sign of det(Hess0) provides us with the parity of the number of eigenvalues of on (0,+oo). 
On the other hand, as we shall see in the next section, the relevant Hessian there is the constrained 
Hessian Even though we do not endeavor to prove it in the present paper, we do expect that the 

parity of the negative signature of A\Xu‘^ of number of eigenvalues of £/ on (0, +oo) coincide, so 
that the results of the current section could be thought of as a direct analogue of the finite-dimensional 
case, with in place of the classical Hessian of the Hamiltonian. The deepest way to prove our 

claim regarding the agreement of those parities consists in examining the Krein signature of eigenvalues. 
Indeed, building on the fact that eigenvectors of are orthogonal for the quadratic form associated with 
A, one may expect to prove that the negative signature of Ayj’^tg is the number of eigenvalues z, with 
Re(z) > 0 and negative Krein signature, and our claim on parity would then follow from the fact that 
eigenvalues with Re(z) > 0 but Im(z) ^ 0 come in pairs. See detailed discussions, precise statements and 
proofs of similar results in [KKS04, BJKll, BJK14] and [KP13, Chapter 7]. 

4 Co-periodic, orbital stability 

4.1 Abstract setting 

We still consider a Hamiltonian system of the form (2), which we relabel here for the reader’s convenience: 

atU = 9ABE.^[U]) , (24) 

with U taking values in B a nonsingular, symmetric N x N matrix, Jtf = ,XJx), and denote 

by ^ the impulse — or momentum — defined by 

^(U) = iu-B-^U. 

Eq. (24) is obviously a system of N (local) conservation laws of order at most three in the spatial 
variable x. Of course we have in mind the more specific forms of B and Jif that are described in Section 
3.1, and correspond to either (qKdV) in the case A = 1, or to (EK) in the case N = 2. In the latter 
case, the first conservation law is of order one, and the second one is of order three as regards the first 
dependent variable, and of order one for the second dependent variable. Section 5 is devoted to a detailed 
investigation of those ‘examples’. Here, we refrain from restricting to any specific form of B and Jif, 
in order to emphasize the crucial ingredients in the proof of co-periodic, orbital stability. The reader is 
referred to Sections 5.1 and 5.2 for an application of our abstract result (Theorem 3 below) to respectively 
(qKdV) and (EK). 

As far as smooth solutions of (24) are concerned, they satisfy at least two additional, local conservation 
laws. One is the conservation of the impulse 4?, Eq. (4), and the other one is the conservation of the 
Hamiltonian Jtf, which explicitly reads 

dtJ^lV] = a^(iEJf’[U] • B EJf’[U] -E ■ B d^EJf[V]) . 

All these local conservation laws have the most important consequence that, along smooth periodic 
solutions to (24), we have 

f Udx = 0, f ^(U)da: = 0, f .^[Uldx = 0, (25) 

dt Jo dt Jo dt Jo 

if 5 denotes the period of those solutions. For a given 5, we call energy space, and denote by Hh a dense 
subspace of (L^(IR/5Z))'^ on which the functional 

[ .^[U]dx 
Jo 
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is (at least) C^. Behind this loose definition, we merely have in mind Hh = for (qKdV), and 

Hh = X L^(IR/SZ) for (EK). Note that the linear functional 

[ IJdx 
Jo 

is automatically C°° on (L^(IR/E!Z))^, by the embedding L^(IR/SZ) ^ L^(IR/SZ), and that the quadratic 
functional _ 

[ i?(U)dx 
Jo 

is also C°° on (L^(IR/E!Z))^ by the Canchy-Schwarz inequality. Therefore, whatever the constants 
(^, A, c) € the functional 


V ^ := [ (^^(U, U^) + c^(U) - A • U + m) dx 

Jo 

is on the energy space Hh. Let us point out that these functionals — in particular ^[U;^, A, c] — 
are invariant under the action of spatial translations U i—>■ U(- + s) on S-periodic functions U, so that 
they are indeed well defined for U viewed as a function on the circle M/SZ. Furthermore, (25) implies 
that ^[•;^, A, c] is preserved along smooth, S-periodic solutions of (24). 

Our main assumptions are the following. 

(HO) There exists an open set O of and a family of periodic traveling profiles U parametrized by 

(/X, A, c) € such that the profile equations in (6)-(7), 

E(Jf’ + ci?)[U] = A , Ljr[U] - c^(U) + A • U = ^, 

hold true, and the mapping (^, A, c) G i-?- (U, S) G C^(K) x K is continuously differentiable, where 
S denotes the period of the profile U. 

(HI) The derivative of the period ^ with respect to the energy level /i, denoted by does not vanish 
on n, and the abbreviated action integral 


0(/x,A,c)= [ (.^(U,U^) + c^(U)-A-U + ^)dx 
Jo 


is such that the matrix 


C := - v^e, 

Va 


is nonsingular for (/x. A, c) G fl, with ^ \ 


(H2) For all in the set of periods achieved on fl, there exists a dense subspace Hh of {LJ 
and an open subset of Hh containing all the profiles U on which the functional 


SZ))^, 


U / M’[\J]dx 
Jo 

is C^, and if we denote by (•, •) the dual product between HL and IHI 2 , there exists a positive number 
a such that 

||U||2^ = (Hess^[U]U,U) + a|jU||i. 

defines an equivalent norm on IHI 2 , uniformly in the parameters defining the S-periodic profile U . 

(H3) For all S in the set of periods achieved on O, there exists a dense subspace Wh of the energy space 
Hh on which the Cauchy problem for (2) is locally well-posed. 


These assumptions are discussed for (qKdV) and (EK) in Section 5. 
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4.2 An index for co-periodic, orbital stability 

Theorem 3. Under the assumptions (H0)-(H1)-(H2)-(H3), for all {p,, A, c) G ri such that 

• the negative signature of C equals the one of the operator A = V\ess{Jtf + c^)[U], 

• the kernel of A is spanned by 

the periodic wave of profile U is conditionally, orbitally stable in the following sense. 

For all e > 0, there exists rj > 0 so that, for all Uq G Wh such that ||Uo — U||hh ^ V> ifT Ihe maximal 
time of existence of the solution U : 1 1 —>■ U(-,t) G Wh to (24) such that U(0) = Uq, then 

inf ||U(-,t)-U(- + s)||H= <e, VfG [0,T). 

In practice, for the cases discussed in Section 5, the fact that the kernel of A is spanned by is a 
consequence of the assumption in (HI) that is nonzero. This is nevertheless a crucial point in the 
proof of Theorem 3, that is why we state it explicitly. Otherwise, the most important and nontrivial 
assumption is 

n(A)-n(C) =0. 

In this respect, the integer n(A) — n(C) may be called an orbital stability index. It is investigated in more 
details in Section 5, where we show in particular its connection with the negative signature of Hess© 
itself, through the remarkable formula 

n(A) - n(C) = n(Hess0) - N . (26) 

Remark 4. At the abstract level of Theorem 3, it is not obvious that the stability criterion n(A) = n(C) 
contains the necessary condition (—1)^ det Hess© > 0 derived in Section 3 (Theorems 1 and 2) for spectral 
stability. However, if we admit (26) for a while, we readily see that a null orbital stability index means that 
n(Hess©) = N, which implies that (—1)'^ det Hess© > 0. For an alternative connection, see Remarks 2 
& 6 . 

Remark 5. In special cases, genuine orbital stability can be inferred from conditional orbital stability. 
This was done^° for example by Bona and Sachs [BS88, Theorem 4] regarding the stability of solitary 
waves in (EKL) with a constant k. Indeed, in this case (EKL) is a semilinear system of PDEs, and a 
bound on the low order derivatives of the energy space yields a bound on higher order derivatives, by 
differentiation of the PDEs and by commutator estimates for the lower order terms. 

Remark 6. Going on with our analogy with the finite-dimensional, ODE case initiated in Remark 2, 
we observe that the above theorem essentially shows that, if the negative signature of is zero (see 

Eq. (27) below) then the wave is nonlinearly stable, while the results of the previous section show that, 
if that negative signature is odd then the wave is spectrally unstable. Even in the Hnite-dimensional 
case, this offers a genuine dichotomy^^ only for planar ODEs. Eor periodic waves, by contrast with 
what happens for the effectively lower-dimensional solitary waves or kinks, it turns out that we never 
have a genuine dichotomy. Nevertheless, by transferring infinite-dimensional conditions on A^^^sg to 
finite-dimensional ones on Hess©, we have come up with the neat following criteria: 

• if n(Hess©) — N is zero then the wave is nonlinearly stable (by Theorem 3 plus (26)); 

• if n(Hess©) — N is odd then the wave is spectrally unstable (by Theorems 1 & 2). 

the sense that it is proved there that solutions starting sufficiently close to the background wave are global in time. 
However the result is still conditional in the sense that Uq is required to have higher regularity — Uq S Wh — than 
afforded by the energy norm || ■ ||h=' 

to considering also the opposite of the ‘natural’ Hamiltonian, if needed. 
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Proof of Theorem 3. It heavily relies on [BGNR13, Theorem 3], in which we proved that 


n(A) — n(A|r^<^) + n(C), (27) 

with 

:= {U e ; /q" U dx = 0, U • Vi?(U) dx = 0} , 

and is a revisited — expanded and more accurate — version of the proof of Corollary 2 in [BGNR13]. 
We proceed in the same spirit as in [DBGRN14, Theorem 7.3]. 

Step 1. Orbital stability within the constraint manifold 

= {UgHh; /o"Udx = /ifUdx, /o"^(U)dx = jf^{lJ)dx}. 

The assumption n(A) = n(C) implies n(A| 7 ’^<g>) = 0 by the formula in (27) recalled above. Since the 
kernel of A is spanned by U^,, this altogether implies the existence of C > 0 so that 

(AV,V) > G||V||i. 


for all V G Hh such that 

/Q"Vdx = 0, /g" V • V^(U) dx = 0, /g"V-U^dx = 0. (28) 

In addition, C is bounded by below by a uniform positive constant when the parameters (A, c) vary in a 
(small) compact subset of A, the projection of 12 onto and /x is implicitly defined as a function of 

(A,c) by the fixed period S (which is made possible by the assumption 7 ^ 0 in (HI)). Therefore, by 
(H2), there exists another positive constant C such that 

(AV,V) > C||V||2^, 

for all V G Us satisfying (28). Indeed, using the equivalent norm on Hh given in (H2), up to augmenting 
a in such a way that 

-cV-B-iV<aV-V, VVgM^, 

(and c possibly varying in a compact set of possible wave velocities), and recalling that A = Hess^[U] + 
cB~^, we see that 

||V||L = (AV,V) - c / V-B~^Vdx + a[ V-Vdx 
Jo Jo 

for V G Hh satisfying (28). 

Now, by Taylor expansion we have 

^[IJ; fi, A, c] - ^[U; /x, A, c] = i (A(U - U), U - U) 

for U G Hh close to U. In the expansion here above, the first order term has vanished because of the 
profile equation E{J^ + c^)[\J] = A, which means that U is a critical point of the functional /x. A, c], 
and the second order term comes from the fact that the operator A = Hess{J(f + c^) is precisely the 
second variational derivative of the functional ^[q /x. A, cj. For those U G Hh close to U that in addition 
belong to we have 

u-u = v + o(|1u-u||hJ, VGru-^, 

which means that V satisfies the constraints in (28) except for the last one. A nowadays well-known trick 
to enforce this constraint is to use translation invariance and the implicit function theorem to prove the 
following. 


< ( 1 +^) (AV,V) 
+ o(||U-U||2j 
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Lemma 1. For all e > 0 we define 

U, = {UgHh; inf ||U-U(- + s)||hh < e} • 

There exists Eq > 0 and a function r : Ugg —>■ K such that for all U G Uej,, 


(U(a; + t(U)) - Vix)) ■ n,(x) dx = 0, |!U(- + t(U)) - U||h, < eq . 


and 

|1U(-+t(U))-U||hh ^0 when inf ||U - U(-+ s)||hh ^ 0. 

As a consequence, for U G Ue;, P| we have by the invariance of /x, A, c] under spatial translations, 
^[U; n, A, c] - ^[U; A, c] = ^[U; A, c] - ^[U; /x. A, c] 

= i(AV,V)dx + o(|!U-U||2j, 

where we have denoted by U the translate U(- + t(U)), so that 

U-U = V + o(||U-U||hJ, VgTu'^, /“V.U,,dx = 0. 

Therefore, we have the lower bound 


^[U;/x,A,c]-^[U;m,A,c] > f ||V||2^ + o(|lU - U||^J 

>f l|U-U||^^, 

for U G Pi 'io and some ei G (0, eq]. 

This is all what we need to use /x. A, c] as a Lyapunov function to show orbital stability within 
Even though this is a classical reasoning, we give it for completeness. For all £ G (0, ei], for all 
U G Uei P such that ||U — U||hh = e, we have 


C 

^[U; /X, A, c] > ^[U; /x. A, c] + — =: to(£) > ^[U; /x, A, c]. 

By continuity of /x. A, c] at U and its invariance under spatial translations, there exists r\ G (0, e] such 
that 

.^[U; /X, A, c] < mis) 

for all U G U^. Therefore, if we take Uq G WsPUr^P*^, and denote by U(t) the solution at time 
t G [0,r) of (24) such that U(0) = Uq, we have 


^[lJ{t);fi, A, c] = ^[Uo; fi, \ c] < m{e) 

for all t G [0,r). Since U(t) belongs to ^ by the conservation of /pUdx and /p=S(U)dx in (25), 
the definition of m(e) and the mean value theorem prevent ||U(t) — U||hh from growing larger than e. 
Indeed, if this happened, there should exist a time t such that ||U(t) — U||hh = £, hence in particular 
U(t) G Uei P^, and therefore 

^[U(t); /X, A, c] > m(e) 

whereas we know that 

^[U(t); /X, A, c] < m(e). 

This proves that U(t) belongs to Ug whenever U(0) G In addition, this rj — as well as the 

Eg, El invoked in the derivation of ry — can be chosen to be the same for parameters (/x, A, c) varying in 
a compact subset of Q on which the period S remains constant, because ^ is uniformly continuous on 
compact sets of Hh x ft, the profile U depends continuously on the parameters (/x, A,c), as well as the 
lower bound C, as already mentioned. This uniformity will be used in a crucial way in the final argument. 
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Step 2. Find a way out of the constraint manifold. 

For fixed (tr, A, c) € let us denote by U the associated profile, of period S, 

U, = {U e Hh ; inf ||U - U(- + s)||h= < e} 

for all e > 0, We still denote by U profiles associated with ‘generic’ parameters {fi, A, c) G There exists 
a neighbordhood of (A, c) in A such that for all (A, c) G B^ and /i = /i(A, c) — prescribed by the fixed 
period S —, the corresponding prohle U belongs to U,,. 

Next, we claim that for (M, P) close enough to (/p“ Uda;, ^(U) dx), say 

|m - /o"ndx|[^ < d(e), |P - J^^(Jl)dx\ < 6{e) 

for some positive d(e), there exists (A, c) G B^ such that the traveling prohle U associated with (^(A, c), A, c) 
satishes 

/o" U dx = M, fg" i?(U) dx = P. 

This follows from the inverse mapping theorem. As a matter of fact, the Jacobian matrix of the mapping 

(A,c)^(/o^Udx,/“^(U)dx) 

turns out to be — C. This is precisely the reason why we have called C the constraint matrix. For more 
details, see [BGNR13, Theorem 3] where we use in a crucial way relations in (10), 

0^ = S, VA0 = -/o“Udx, 0e = /o"i?(U)dx. 

Since by (HI) C is assumed to be nonsingular, the inverse mapping theorem does apply and this proves 
our claim. 

Conclusion. We now have all the ingredients to complete the proof of Theorem 3. Let us take e G (0, ei], 
where ei is introduced as in Step 1 for (y. A, c) and its neighbors. This is where we use uniformity: by 
Step 1, there exists rj G (0,£] so that, for all (A, c) G B^^, if Uq G Wh is such that 

inf II Uo - U(- + s)||He <V, Jo Uo dx = /g" Udx , jf ^(Uq) dx = /“ ^(U) dx , (29) 

with U being the traveling profile associated with (/x(A, c). A, c), then U(t) G Ue/2 for all t G [0, T), where 
T is the maximal time of existence of the solution U(t) of (24) such that U(0) = Uq. 

By continuity of the mapping 

U G Hs ^ Uo U dx, Jf i?(U) dx) 
there exists ( G (0,77/2] such that ||Uo — U||hh < C implies 

|/“Uodx - /o"ndx||^^ < <5(77/2), |/“^(Uo)dx - /o"i2(n)dx| < 6^/2), 

where <5 is the function involved in Step 2. This implies by Step 2 the existence of (A, c) G B^/2, its 
associated profile being U G 11^/2 ®'^ch that 

Jo U dx = /p" Uo dx, /p" ^{V) dx = J^ ^(Uq) dx . 

Therefore, if Uq G Wh is such that ||Uo — U||h 3 < C, there exists a profile U such that we have (29) 
— by using the triangle inequality to achieve the first condition. By Step 1 this implies that U (t) G Ue/2 
for all t G [0,r), and thus U(t) G by the triangle inequality again. This proves the orbital stability of 
the traveling wave (x, t) <-)■ U(x — ct). □ 
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Before concentrating separately on the cases N = 1 and N = 2, let us point out a general identity on 
the constraint matrix and the Hessian of the abbreviated action integral. 


Proposition 1. Let 0 : H —>■ M 6e a function of (/x, A,c) € Lt, an open subset o/M x x M, 
such that the second order derivative does not vanish, and define the continuous function C : H —)■ 

^(N+l)x{N+l) Jjy 


c = 


V0^(g) V0^ 

©MM 




where V stands for the partial gradient with respect to all but the first independent variable p. Then there 
exists a continuous mapping P : H —> SL7v+2(IR) such that 


P (Hess0) 



In particular, we have 

det(Hess0) = (—0^^ det C . (30) 

Proof. The matrix relation is a matter of elementary operations on rows and columns, which give the 
result with 


P = 


Eq. (30) is a straightforward consequence of that relation. 


1 

0 

e^V0^ 

In+1 


□ 


As a consequence, we see that n(Hess0) = n(—C) if 0^^ > 0 and n(Hess0) = n(—C) + 1 if 0^^ < 0. 
This is the key to the proof of the identity in (26), together with a Sturm-Liouville argument as we 
explain in the next section. 


5 Examples 

5.1 Quasilinear KdV 

We use here the notational convention introduced before Theorem 2, namely, for a given periodic wave 
solution to (qKdV) of spatial period T and profile v, 


9{p,, A, c) = J {£{v, v„,) + cr^v) — Xv + p)dx , 


/o 

a = Hess(.e+ cxj)[v\ , 

Accordingly, we denote by I the Lagrangian such that 

rX 


{v) = \v\ 


0{p,,X,c)= / {£[v; X, c] + fj,) dx, 
Jo 


'0 

and whose second variational derivative is a, that is 


£[v; A, c] ;= £{v, v^) + \cv‘^ — Xv. 

The assumption in (HO) regarding the existence and parametrization of periodic wave profiles is easily 
met when the energy is of the same form as in (3), 

£iv,Vx) = f{v) + Ik{v)vI , 
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with smooth functions / and k such that k{v) > 0. Indeed, the profile equation Lf [u; A, c] = /r then reads 


- f{v) - + Xv = fi. 

Remarkably enough, the associated phase portrait in the plane {(f,!!); v := Vx\/k{v)} does not depend 
on K, and consists of the level sets 


{few); -/fe) - + Xv = fi} ■ 

When /i is varied, these level sets exhibit saddle points {v, 0) where the potential 

W{v; A, c) := —f{v) — + Xv 

achieves a local maximum, and center points (v, 0) where 1^(-; A, c)) has a local minimum. In other words, 
since /' = —p, {v, 0) is a saddle point if p{v) — cu +A = 0, p'{v) < c, and a center point if p(v) — cv + X = 0, 
p'{v) > c. In particular, if p is convex, a family of periodic wave profiles is found as soon as we have a 
pair made of a saddle point (us, 0 ) and a center point (uoiO) such that there exists v® with Vg < vq < 
and 

W{vs]X,c) = 1f(v’^\X,c ), 

which amounts to an equal area rule on the graph of p (that is, the areas in between the graph of p and 
that of the affine function u i—>■ cu — A, for u € [ws, and for v € [wqj are equal; this can be observed 
on Figures 1 and 2 hereafter). 

By the implicit function theorem, the roots v of p(v) — cv + X are smoothly parametrized by (A, c) as 
long asp' — c does not vanish, which means that saddle points and center points are smoothly parametrized 
by (A, c), and the zeroes of W (•; X,c) — p are smoothly parametrized by (/i. A, c) away from critical points 
Vs{X,c) and uo(A,c). By smoothness of the flow of ODEs, this shows that periodic wave orbits found 
inside homoclinic loops are smoothly parametrized by {p, A, c). 

Let us give a more precise situation in which (HO) is satisfied. It is chosen in order to include many 
of the examples we have in mind, and in particular power laws p{v) = v'^, 7 > 1 , and p(v) = v~'^, 7 > 0 . 
(The more complicated van der Waals law has been investigated in earlier work, see [BGNR14].) Let us 
point out in passing that, as far as profile equations are concerned, the special case p(v) = which 
corresponds to the ‘standard’ KdV equation, and p(v) = v~‘^, which corresponds to a shallow-water type 
of pressure law^^, are closely related. Indeed, the profile equation for p{v) = readily amounts to a cubic 
potential W, while the profile equation for p{v) = v~‘^ also amounts to a cubic potential after multiplying 
it by V — and modifying k accordingly (one should however pay attention to the fact that this operation 
alters the status of the parameters (/x, A,c), and in particular that of p, which becomes like a Lagrange 
multiplier instead of being an energy level). Even without this trick, the phase portraits look similar, see 
Figures 1 and 2 to compare the two situations. 


^^See §5.2.3 for an explanation. 
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Figure 1: Quadratic nonlinearity, associated potential and phase portrait. 
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Figure 2: Shallow-water type nonlinearity, associated potential and phase portrait. 
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Proposition 2. Assume that k : / — )■ (0, +oo) is C^, and that p = —f : / — > M is with p" > 0 on the 
open interval / C M. We denote by V : J := p'{I) —>■ I the inverse mapping of p', and 

I : c G J 1 -^ cV{c) — p{V (c)), 

I- : c G J i-G lim (cv — p(v)), 

i7\iirLf I 

1+ : c G J i-G lim (cv — p(v)), 

V /^sup I 


and eonsider the open sets 

r := {c G J; l-(c) < c, l+(c) < 1(c)} , 

A := {(A, c) G M X r ; max(/_(c), l+(c)) < A < ^(c)} . 
Then there exist mappings 


Vs '■ A ^ {v G IV <V (c)} and oq : A —>• {n G /; v > V(c)} 
such that for all (p, A, c) G fl, 

n := {(p, X, c) G R X R X r ; (\,c) G A, p G (W (no(A, c); A, c), W (vs(X, c); A, c)))} , 
A:={(A,c)gA; lim W(v; X, c) > W(vs(X, c); X,c)} , 

V /^sup I 

there is a unique solution v to 

Ik(v)vI - f(v) - ^cv^ + Xv = p, vjO) = 0, 

that is periodie, and it is in parameters (p, A, c). 

The proof is based on the remarks made above and some elementary analysis of the variations of 
W(-]X,c) for (A, c) G A, see Table 1 below. Details are left to the reader. 


V 

V\ 

Vs 

V2 

Vo 

V3 v'^ 

w 

+ 

0 

- 

0 

+ 

w 

T 


\ 

T 

\ 

Po 

T 


Table 1: Variations of potential 

Proposition 2 applies in particular to 

• p('v) = w'*', 7 > 1, / = (0,+oo), J = (— 00 , 0 ), V(c) = 

l-(c) = 0 , l+(c) = —CO , 1(c) = ( 1 / 7 ^/^^“^^ — , 

• , 7 > 0, / = (0, + 00 ), J = (0, + 00 ), V(c) = (—^/c)^Ai+A^ 

l-(c) = —00 , l+(c) = —00 , 1(c) = (yl/lv+l) _ J^y'^7/(7+l)^('_gp/(7+l) _ 
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This is what we may say on (HO). Regarding (HI), it turns out to be equivalent to 

det(Hessd) 7 ^ 0. 

Indeed, 7 ^ 0 is exactly the first condition in (HI), and as soon as it is satisfied, the equivalence 
between detc 7 ^ 0 and det(Hess0) 7 ^ 0 readily follows from Proposition 1, in which Eq. (30) reads, with 
our present notation (and = 1 ), 

det(Hessd) = det c. 

We may also derive this relation by writing explicitly 

\ 


and by computing 




( 

0\\ 



Q\c 


1 





0\IJ. 


Ofifj. 

1 

dc\ 



dec 




dcn 



Ocfj, 

d\\ 



S\c 

0fj,c 




0,, 


d\ll 



= 

0c\ 



Occ 

Ofic 



dcp. 

Ofifj. 



0^^ 
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det(Hessd). 


The investigation of whether 0^^, the derivative of the period with respect to the value of the ODE 
Hamiltonian vanishes is a classical topic in Hamiltonian dynamics, see for instance the recent paper 
[GV14] and references therein. A criterion ensuring that this derivative is positive was given in particular 
by Chicone [Chi87]. It merely reads {W/> 0, for a Hamiltonian of the form ^v‘^ + W{v). Despite 
its simple form, it is not easy to check analytically. We have chosen to rely on numerical experiments in 
Section 6 to rule out the critical cases in which 9^^ would be zero. Regarding the zeroes of the determinant 
of Hessd, far from distinguished limits — see [BGMR15] for a discussion of those —, we are not aware of 
any general analytical result. This is also investigated numerically in Section 6 . 

Concerning assumptions in (H2), we claim — recall that the period of waves is denoted by T instead 
of S here — that the space H-f = i7^(IR/TZ) is a convenient choice as soon as 


^{v,Vx) = f{v) + hK{v)vl 


with f '■ I 


and K : I —>■ {0, + 00 ) of class on an open interval I. Indeed, in this case, the mapping 

rT 


E -.v ^ 


£{v,Vx) dx 


is well-defined on the open subset of i7^(IR/TZ) made of v with values in J — thanks to the embedding 
iJ^(IR/TZ) ^ — and twice differentiable with 

T 

dE{v) ■ h = f {f{v)h+^K'{v)hvl + K{v)vxhx)dx,yhGH^{R/TZ), 

Jo 

d^E{v) ■ {h, k)= [ {f{v)hk + ^K"{v)hkvl + K'{v){hkx + kK) + n{v)hxkx)dx ,yh,ke H^{R/TZ ), 
Jo 

for all V G E^(R/TZ) with image in the domain of definition of / and k. If in addition v G C^, as is the 
case for traveling profiles, we may integrate by parts in the formulas above, and recognize the variational 
derivatives of £. As a matter of fact, if v G we find that 

dE{v)-h= [ hE4v]dx,yhGH^{R/TZ), 

Jo 
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d?E{v) ■ {h, k) = (Hess^H/i, k) ,\/h,k& iJi(K/TZ), 

where 

E^[w] ;= dy£- d^{dy^£) = f'{v) + Ik{v)vI - dxiKiv)va;) G C° , 

and 

Hess£[t;]/i := {d'^£- d^{dyy^£)) h - d^{hdl^£) = (/"(w) + \k''{ v)vl - d^{K\v)v^)) h - d^{K{v)h^) 

belongs to iJ“^(IR/TZ) = (i/^(IR/TZ))' for all h G i/^(M/TZ), and (•,•) denotes the dual product. In 
particular, Hess.e[v] is a Sturm-Liouville operator, of the form —d^Kd^ + <?, with T-periodic coefficients, 

K = k{v) ec^, 0 < Ko< K < Ki, 

q = /"{v) + ^k"{v)vI - da;iK{v)Vx) G C° , ||(7|U“ < ^0 • 

Therefore, for all h G i7^(M/TZ), we have 

{Hess£[v]h, h)a\\h \\'^2 = ( {Kh^{qa)h'^)dx , 

Jo 

and if we choose for instance a = Oq + 1; we have 

l-r i-r 

min(iGo 5 1) / +/i^)da; < (Hess£[u]/i, h) + a||/i |||2 < max(iGi, 2 q;) / +/i^)da;. 

Jo Jo 

This shows the equivalence of norms requested in (H2). 

Finally, the main assumption in (H3) is satisfied at least in i7^(IR/TZ) for s > 7/2. Indeed, the 
following theorem is proved in a forthcoming paper [MieI5]. 

Theorem 4. Let s be an integer such that s > 112. If p = —/' : I ^ R is and k : I —)■ (0, +oo) is 
then for all T > 0, vq G i7'*(M/TZ), the image of vq being in I, there exists T > 0 and a unique 
V G C([0, T); i7^(IR/TZ)) solution to (qKdV) with e = f{v) + ^k{v)v'^, and vq v is continuous. 

In the special case when k is constant, well-posedness is known to hold true with much lower regularity, 
for example s > —1/2 is sufficient for the classical KdV [KPV96]. However low-regularity results rely on 
dipsersive effects to control nonlinear terms and are thus strongly model-dependent. This is notoriously 
still a held of intense research^^. 

Theorem 5. Under the assumptions of Theorem 4, if (HO) is satisfied (in particular if the assumptions 
of — and thus — Proposition 2 hold true), and if, for some (/i. A, c) G O, 

(s) det(Hess0) 7 ^ 0, n(Hess0) = I, 

is satisfied, then the periodic wave associated with {p, A, c) is conditionally, orbitally stable in i7^(IR/TZ). 

Proof. The assumptions (H2)-(H3) are implied by those of Theorem 4, while, as explained before, (HI) 
is equivalent to 

det(Hess0) ^ 0. 

So our main assumptions (H0)-(H1)-(H2)-(H3) are met here. 

Furthermore, as announced before, the assumption that the kernel of a is spanned by v„. is automat¬ 
ically satished. In fact, Lemma 2 stated below applies to the present operator a, by using the potential 
W = W{v, A, c) introduced at the beginning of this section. We thus infer at the same time that the 
kernel of a in L^(IR/TZ) is spanned by v„,, and how to compute the signature of a. Therefore, to apply 

^^Incidentally we point to the attention of the reader the attempt of keeping track of latest known results — including 
local-well posedness proofs — for various dispersive equations on http://wiki.math.toronto.edu/DispersiveWiki/ 
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Theorem 3 and thus complete the proof of Theorem 5 it suffices to check that our assumptions imply 
n(a) == n(c). 

To do so, we can actually bypass the computation of n(c) by using Proposition 1. With our current 
notations, this algebraic proposition shows indeed that 

n(Hessd) = n(—c) if 9^^ > 0, n(Hess0) = n(—c) + 1 if < 0, 

while we know from Lemma 2 that 

n(a) = 1 if > 0 , n(a) = 2 if < 0 . 

Therefore, regardless of the sign of — as long as this number is nonzero —, we have 

n(Hess6*) = n(—c) + n(a) — 1. 

Furthermore, since the 2x2 matrix c is nonsingular — as already justified —, 

n(-c) = 2 - n(c), 

so that the previous formula equivalently reads 

n(Hess0) = 1 + n(a) — n(c). 

This is the announced identity in (26) in the special case = 1. 

As a consequence, the stability index n(a) — n(c) vanishes if and only if 

n(Hessd) = 1. 

In this respect, the set of conditions in (s) are the optimal ones enabling us to apply Theorem 3 to 
(qKdV). □ 

Lemma 2. Assume that k : I —> (0, +oo) and W ■. I are on some open interval I and such that 
the Euler-Lagrange equation E'^[v] = 0 associated with the energy 

^ : {v,Vx) i-t ^k(v)v‘^ + W{v) 

admits a family of periodic solutions v taking values in I, parametrized by the energy level p, = L^[u], for 
p & J, another open interval. If we denote by T the period of v, and assume that T^, its derivative with 
respect to p, does not vanish, then the self-adjoint differential operator a := Hess^ju] has the following 
properties: 

• the kernel of a on L^(IR/TZ) is the line spanned by p„.; 

• the negative signature n(a) of a is given by the following rule: 

* > 0 then n(a) = 1, 

* < 0 then n(a) = 2. 

Since 

a = -d,^K{y)dx + {v)v „,), 

is a Sturm-Liouville operator with periodic coefficients, the main argument in the proof of Lemma 2 
relies on Sturm’s oscillation theorem, as for instance in Lemma 1 in [BJKll], which concerns the case of 
a constant k. The detailed proof is postponed to Appendix B. 

Note that since Hess0 is a 3 x 3 matrix, its determinant cannot be positive if n(Hess6*) = 1. So we 
could equivalently replace the second condition in (s) by det(Hessd) < 0. 

If the first two conditions in (s) are readily amenable to computations, the verification of the third 
one, n(Hess6*) = 1 demands some slightly more sophisticated algebraic work. It can be interesting to have 
more explicit conditions, in particular to compare with earlier results. If we assume moreover that 2x2 
determinant 9\\9^^ — does not vanish, we can see by an elementary count of sign changes in the 
principal minors of Hessd that (s) stems from having either one of the following sets of conditions 
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(sl)0^^>O, ^0, det(Hess0) <0, 

^A/i 

(s2)0^^<O, <0, det(Hess0) <0. 

PA/i t/fifi 

For convenience, the reader may refer to Table 5 in Appendix A, in which (si) corresponds to the second 
and third rows, and (s2) to the 7th. From the same table we see that, in terms of the constraint matrix c, 
(si) corresponds to n(c) = 1 , and (s2) to n(c) = 2. Note however that our conditions in (s) are slightly 
more general than the prescription of ((si) or (s2)), in that they do not require — 9\^ ^ 0 . 

In the special case when k is constant. Theorem 5 was essentially already known in the case ^aa^/x/x ~ 
7 ^ 0, and proved in a slightly different manner in [BJKll]. Indeed, orbital stability with respect to 
co-periodic perturbations is essentially a consequence of [BJKll, Theorem 1], under the assumption 

6 '/xm 7 ^ 0 , 7 ^ 0 , det(Hess 6 ») 7 ^ 0 , n(Hess 6 ») = 1, 

which is equivalent to ((si) or (s2)). An earlier, similar result was shown by Johnson [Joh09], under the 
more restrictive assumption 

<0, det(Hess0) < 0. (31) 

We should mention that these results by Bronski et al dealing with (gKdV) — and not its quasilinear 
version (qKdV) — yield genuine orbital stability, and not only conditional orbital stability. This is 
because the Cauchy problem for (gKdV) is much better understood that for (qKdV) with a nonconstant 

K. 

Remarkably enough. Table 5 in Appendix A shows that the matrix Hess0 cannot be nonnegative. 
Indeed, by our identity on negative signatures (see Eq. (27)), we must have n(a) > n(c), so that the — 
purely algebraic — situation in the first row of Table 5 cannot occur. 

5.2 Euler—Korteweg 

5.2.1 Eulerian coordinates vs mass Lagrangian coordinates 

Before investigating how to apply Theorem 3 to the Euler-Korteweg system, let us come back to the 
equivalence between its formulation in Eulerian coordinates (EKE), and its formulation in mass La¬ 
grangian coordinates (EKL). This equivalence works as long as we deal with states away from vacuum, 
and more precisely, with densities p that are bounded and bounded by below by some positive constant. 
It is based on the fact that the continuity equation 

dtp + dx{pu) = 0 (32) 

is equivalent — for a: € M, and t in some interval too — to the existence of a function Y = Y{x,t) such 
that dxY = p and dfY = —pu. Denoting by (y = Y(x,t),s = i) the new coordinates obtained this way, 
and introducing 

v{y, s) = l/p{x, t ), w{y, s) = u{x, t ), 

we see that (EKE) and (EKL) are equivalent provided that (S{p,px) = p£{v,Vy). (Here we focus on 
smooth solutions, but this equivalence is also known to hold true for weak solutions when S’ depends only 
on p, that is, for the usual Euler equations.) The change of coordinates {x, t, p, u) 1 —)■ {y, s, v, w) is clearly 
nonlinear, and also nonlocal since the new ‘independent variable’ y actually depends on the integration 
of the dependent variables p and pu. Nevertheless, there is also an equivalence between travelling wave 
solutions of (EKE) and (EKL). As shown in [BG13], {p,u) = {p,u){x — at) is a travelling wave solution 
to (EKE) if and only if {v,w) = {y,w){y + jt) is a travelling wave solution to (EKE), along with 

p{u-a)=j, p{x) = l/v(Y{x))), u{x) = w{Y{x)), Y'{x)=p(x), 
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or, equivalently, 


w-jv = a, v{y) = l/g{X{y))), w{y) = u{X{y)), X^{y)=y{y). 


In particular, if {g , m) is ^-periodic then (u, w) is T-periodic with 


T = / g{x)dx, - = / v{y)<^y- 

Jo Jo 

Furthermore, we can see that the abbreviated action integral does not depend on the chosen formulation. 
Before checking this, let us perform preliminary work on profile equations. The system (EKE) fits our 
abstract framework with 

U = ( ^ ) , J^[U] = S{g, g,) + \gu^ , ^(U) = -gu , 
so that the profile equations in (6)-(7) can be written as 


{dp^){g,gx) - {{dp^(S’){g,gx))x + - au = -A, 

gu-ag = j , (33) 

Bx idp^(g’){g,gx) - ^{g,gx) - ^gy? + a gu - Xg + ju = g. 


The notational choice for the Lagrange multiplier in the right-hand side of the second equation here above 
is dictated by the relation we already have in mind, and we denote by —A the first Lagrange multiplier 
for convenience, because it is going to be identified with the energy level in the EKL profile equations. As 
a matter of fact, it results from Theorem 1 in [BG13] that the profile equations in Eulerian coordinates, 
as written above, are equivalent to 


{ {dy£){v,Vy) - {{dy^£){v,Vy))y - jW = -/T , 

w-jy = a, (34) 

Vy {dy^£){y,yy) - e{y,Vy) - Jyf' + jyw - gy + aw = A. 

We recognize here the profile equations for (EKL), which fits our abstract framework with 


U = ^ J , = e{v, Vy) + , ^(U) = vw . 

So there is an almost perfect symmetry in these profile equations, where we can see that j and cr exchange 
their roles, as well as A and g, when we go from (EKE) to (EKL) or vice versa. This is why, even though 
this might seem confusing at first glance, we avoid introducing any additional piece of notation. The fact 
that we have a single abbreviated action integral for both (EKE) and (EKL) is now clear because the 
change of variables is such that dy = gdx, or equivalently dx = ydy, so that 


'■= / {S{B,gx) + \bu^ - (ygu +Xg - JU +g)dx (35) 

Jo 

T 

(.e(u, Vy) + — jyw_ + gy_— aw -I- A) dy. 

The third of profile equations in (33) and (34) show that 0 is indeed the abbreviated action for both 
of them. In addition, we have the following correspondence between the ‘abstract parameters’ {g,\,c), 
as introduced in profile equations (6)-(7) and used in the abbreviated action (9), and the ‘practical 
parameters’ {g,X,j,a) for (EKE) and (EKL). 



abstract 

y 

Ai 

^2 

c 

EKE 

y 

-A 

J 

a 

EKL 

X 

-y 

a 

-j 


Table 2: Notation for wave parameters 
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Remark 7. On the one hand, the occurrence of minus signs in Table 2 is not a real issue for applying 
our theory to (EKE) and (EKL), because both the quadratic form associated with HessO and the one 
associated with the constraint matrix C are invariant under the symmetry v —v in any dependent 
variable v of 0, and these quadratic forms are, together with the derivatives and Ta, the only objects 
which govern our stability criteria. On the other hand, the fact that the roles of parameters change when 
we go from (EKL) to (EKE) will have to be addressed carefully. 

The invariance of spectral stability properties when going from Eulerian coordinates to mass La- 
grangian coordinates, even though it seems very natural, is not obvious either as the original conjugacy 
occurs through the change of coordinates {x, t, p, u) i—)■ (y, s, v, w) that is nonlinear and nonlocal. Nonethe¬ 
less, there is a kind of ‘conjugacy’ between systems that are obtained by linearizing (EKE) and (EKL), in 
moving frames, about {p, u) and (w, w) respectively. This in turn enables us to prove that the existence 
of an unstable mode in either one of these systems implies so for the other. In order to point out this 
conjugacy, we first reformulate (EKE) and (EKL) in moving frames associated with the waves. Let a 
and j be fixed, and consider the new dependent variables 

q := p{u-a) - j , z := w - jv - a , 
and the new independent variables 

^:=x-at, (:=y + js. 


Then the system (EKE) is equivalent to 


(EKEj) 


dtp + dxq = 0, 

dtq + do,<i> = 0, $ = <^{p,p^,pa;x,q;j) := (q + /p + pEt^[p] + LS‘[p ], 


where dt now denotes the partial derivative at constant and for convenience we have substituted again 
X for The second equation in (EKEj) here above comes from the conservation law for the impulse, as 
in (4), which reduces here to 


dt{pu) + dx{pv? p E(o’[p] -|- L<l’[p]) — 0 , 


from which we have subtracted a times the continuity equation (32). Similarly, the system (EKL) is 
found to be equivalent to 


(EKLj) 


dsV = dyZ, 

dsZ = dy'ii, 'i' = 'ii{v,Vy,Vyy,z;j) :=E4'^]-2jz-fv. 


Here above, dg stands for the partial derivative at constant and for convenience we have substituted 
again y for (. Of course, as is the case for (EKE) and (EKL), systems (EKEj) and (EKLj) are equivalent, 
as long as smooth solutions with positive and bounded densities p and volumes v are concerned and 
provided that £‘{p,px) = pe(y,Vy). More precisely, the change of coordinates (x,t,p,q) i—>■ {y,s,v,w) is 
given by 

dy = pdx-qdt, s = t, v{y,s) = l/p{x,t), z{y,s) = q{x,t)/p{x,t), 
or equivalently by 

dx = udyzds , t = s, p{x,t) = l/v{y,s), q{x,t) = z{y,s)/v{y,s). 

By construction of (EKEj) and (EKLj), the travelling wave solutions to (EKE) and (EKL) considered 
above become stationary solutions to (EKEj) and (EKLj), which read (p,0) and (w, 0) respectively. We 
can now show the following. 
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Theorem 6. Assume that g and v are smooth functions on M, bounded and bounded by below by positive 
constants, related by 

g{x) = l/v{Y{x))) = Y'{x), 

and such that {g,0) and (v, 0) are stationary solutions to (EKEj) and (EKLj) respectively, for some real 
number j. Then by linearizing (EKEj) and (EKLj) about {g,0) and (v, 0) respectively, we receive systems 
whose spectra are identical. 

Proof. Let us simply call (ffij) and {TLj) the linearized systems of (EKEj) and (EKLj) about (p,0) and 
(v, 0) respectively, and denote them in abstract form as 

m) 

(fLi) 3.())=^.(")^ 

Our aim is to show that the differential operators and Afr are isospectral. By translation invariance of 
(fEj) and (^Lj), we already know that :/ = 0 is an eigenvalue of both and .2^, associated eigenvectors 
being and (wy,0)^ respectively. From now on, we take a nonzero complex number v, and aim 

at showing that it belongs to the spectrum of Afs if and only if it belongs to the spectrum of 2fL- By 
spectrum of the operator .2 e, which is a differential operator in x with S—periodic coefficients, we mean 
the whole spectrum in the space of square integrable functions, which is known to be the collection of 
complex numbers e such that there is a nontrivial (p, q) satisfying 

= iz , 3aeR, {p,q){x + E)=e’-°‘{p,q){x),yx eR. (36) 

Here above, a is called a Floquet exponent. Note that the spectrum of in the space of square 
integrable S-periodic functions corresponds to those v for which 0 = 0. This is the case for i> = 0, 
since (pa:,0) is S-periodic. Of course, the spectrum of enjoys a similar characterization, which is the 
existence of a nontrivial {v, z) such that 


, 3aeR, {v,z){y + T)=e^^{v,z){y),\fyeR. (37) 

The idea is to show a one-to-one correspondence between nontrivial (p, q) satisfying (36) and nontrivial 
(u,i) satisfying (37), for the very same values of (jz, a) with e 0. This can be done by first returning 
to nonlinear systems. 

Solving the original, nonlinear system (EKEj) for some perturbation of p that is parametrized by say 
e as initial data, we receive a family of solutions (p, (?) = {p,q){x,t;e) of (EKEj) parametrized by e such 
that 

{p,q){x,t;0) = (p(a;),0), 

and that {p,q) : {x,t) >->■ {p^,q^){x,t-,Q) solves (ffij). Furthermore, introducing Y = Y{x,t-,e) such that 
Y{x,t;0) = Yfx) and d^Y = p, dtY = —q, we have 

p{x,t;e) = l/v(Y(x,t;e),t ;€), q{x,t) = z{Y{x,t;e),t-,€)/v{Y{x,t;e),t;e) (38) 

where (w, z) = (v, z)(y, s; to) is a family of solutions to (EKLj) parametrized by w. (We use here different 
notations for parameters e and lo for the same reason as for times t and s, that is, in order to avoid 
confusion about partial derivatives.) This implies that 

(i’,2)(p,s;0) = (p(j/),0), 
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and that {v,q) : {y,s) i-)- {vu:,qu:){y,s;0) solves (t'Lj). Assume moreover that the family (p, ( 7 )(-, 0; e) is 
chosen such that (po,qo) = (Pe, 0; 0) satisfies (36). Then 

(p,q) = e‘'*(po,qo), 

and we claim that, similarly, 

(v,z) = e‘"^(vo,io), 

with (wo,io) = satisfying (37). 

In order to prove that claim, let us first note that, by the chain rule applied to (38), 

p(x, t) = -{v{Y_{x),t) + Y (x, t) Vy{Y_{x)))/y{Y_{x)f , q{x, t) = z(Y{x), t)/v{Y{x)), (39) 

where Y(x, t) := Y^{x, t; 0). Now, by differentiating d^Y = p, dtY = —q, we get that Y^ = P,Yt = —q. By 
the first row in (iEj) and the fact that (p, q) depends on t as a linear function of e"^*, we have vp = —qx, 
and therefore, we find that Y = —q/v. This relation is the key to the claimed conjugacy, because it 
enables us to rewrite (39) as 

s) = -y(yfp(X(y), s) + q(X(y), s) Vy{y)/v , z{y, s) = y{y)q{X{y), s) , (40) 

with A = Equation (40) obviously implies that 

v{y, s) = e''^ i-y{y)po{X{y)) + 9 o(A(j/)) yy{y)/v) = e‘'%o{y), 
z(y, s) = e"'"’{yiy)qoiX{y)) = e''^zo{y). 

In addition, we observe that v, z given by (40) in terms of u. A, p, q, and v, are bounded functions of 
p if p, <7 are bounded functions of a; — because y and its derivatives are bounded. Furthermore, v, z 
are T-periodic in y if p, q are S-periodic in a: — because A(p + T) — A(p) = S, and more generally, if 
(p, q){x + S, t) = e*“(p, q)ix, t) for all x, t, then (v, z){y + T, s) = z){y, s) for all y, s. 

Similar computations can be performed in the other way round. Indeed, we can find (y, z) solving 
(EKLj) from {p,q) solving (EKEj) through the following formula, analogous to (39), 

'viy,s) = -{KKiy),s) + x{y,s) Px{X{y)))/p{X{y)f ,z{y,s) = q{X{y),s)/p{X{y)) , (41) 

and the key relation is A = zjv. Using in addition that p = 1/v, Px != —Vyly, we thus see that, as 
expected, (41) is equivalent to (40). □ 

Remark 8. Theorem 6 shows in particular that the operators just obtained by linearizing (EKE) and 
(EKL) in moving frames, but in the ‘original’ dependent variables (p, w) and {v,w), 

■sYe = —Da; [ L, , .£ I and sYl —Ytyi^, K . , 

Hess<^[p] jIp ) ^ y Hess£[p] -j J ’ 

are isospectral. Furthermore, we can infer from its proof the relationship between the eigenfunctions of 
sYe and jYl associated with nonzero eigenvalues. Indeed, recalling that 

q := p{u-a) - j , z := w - jv - a , 

we have 

9 = {j/e)p + eu, z = w- jv, 

which yields, by substitution in (40), 

'V 

V = y^{—l + jV}p + Lu , w = j^'i^LpY {1+ jV}u, L:= — 

vy 

(where we have omitted to write the independent variables for simplicity). The practical ‘conjugacy’ 
between sYe and sYe is thus far from being trivial. 
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Another natural question is the relationship between the Hessians of the constrained energies 

~i) ■ 

Interestingly, both of these matrix-valued operators are linked in a rather simple manner to scalar oper¬ 
ators. Indeed, using notation as above, 

4= {j/e)p + eu, z = w-jv, 

we see that 

( h ) ■ ^^ ( h ) = + q^/p , a.E := Hess«?[^] - f , 

• Al ( ^ ) = v&lv + , ar ■= Hess£[r;] - f . 

These expressions have the striking consequence in terms of negative signatures that 

n(As) = n(aB), n(Ai) = n(ai). (42) 

Both aE and a^ are scalar Sturm-Liouville operators with periodic coefficients. We thus have a simple 
criterion to compute their negative signature (see Lemma 2 stated in Section 5.1). 

Regarding the relationship between the scalar operators aE and a^, one may check that 

g aEP -b K - p^^ p) = - aLV , K := , 

or equivalently 

vaLV + K {Vy Vy - Vyy v) = - aEp , K '.= 9^^, 

for 

v{Y{x)) = -p{x)lp{xf , 

which amounts to substituting 0 for X in (41). These relations are — fortunately — consistent with 
the fact that aEPx and aEVy vanish simultaneously, but they do not show a relationship between the 
signatures of a^; and a^. This is actually no surprise because we expect that it is the stability indices 

li; := n(A£;) - n(C£;) = n(a£;) - n(C£;) and II := n(AL) - n(CL) = n(aL) - n(CL) 

which vanish simultaneously (see Remark 9 here after for a variational point of view on this question), 
and the negative signatures of the constraint matrices Ce and Cl have no a priori reason to coincide. 
Indeed, recalling from (35) the definition of 0(/i, A, j, cr), and interpreting the abstract definition of the 
constraint matrix C in (HI) with the present notation (see Remark 7), we find that 

®AA ®/iA ®Aj ®A(7 

©A/j, ®A/j, ^A/i. 

BjA B/iA Bjj Bjcr B^(j 

Bjp B^^ Bj^ B^p Bj^ B^p 

Bcta B^a Bctj B^j Bo-(7 B^ct 

Bcr/i B^^ B^/j, B^^ B(j/j, B^^ 

B^^ B^^ B^cr Ba(t B/j,j Baj 

B^iA Baa B^a Baa B^a Baa 

Bo-/i, Ba// B(J(7 BA(T BcTJ Baj 

Bcta Baa B^a Baa B^-a Baa 

Bj^ Ba/j, Bjct Bact Bjj Baj 

Bja Baa B^a Baa B^a Baa 
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Obviously, the matrices Ce and do not involve exactly the same minors of HessO, and may therefore 
have different negative signatures — according to our algebraic computations in Appendix A, it might 
happen that n(CL) = 2 and n(C£;) = 0. 

Remark 9. From (35) we see that 

0(M,A,j,cr) = A, j,cr] = A, j,CT] 


^E[p,u]^x,X,j,a] := / {S{p,p^)+ \pu'^ - apu +Xp-ju +p)dx , 

^0 

= / {e{v,Vy)-\-^'uP'— jvw-\-liv — aw-\-\)dy . 

Jo 

The nullity of the stability index 1^; is a necessary condition for ^e['': P,X,j,a] to have a local minimum 
at {p , u) under the constraints 

dx = E, p{x)dx = fg"" p(x)dx , 

^ ^ = .. ( 45 ) 

/q“ u(x)dx = /“ u(x)dx , fg^(pu)(x)dx = fi^“(pu)(x)dx . 

Similarly, the nullity of the stability index 1^ is a necessary condition for P, X, j,a] to have a local 
minimum at {v, w) under the constraints 

lo dy = T, v{y)dy = y{y)dy , 

fo w(y)dy = fg w(y)dy, (vw)(y)dy = (yw)(y)dy. 

We have inserted the apparently trivial, first condition in (45) and (46) in order to point out that these 
sets of constraints are actually equivalent under the change of variables 

p(x) = l/v(Y(x ))), u(x) = w(Y(x )), Y'(x) = p(x ), 

which happens to also ensure that 

^Llv,w;p,X,j,a] = ^Elp,u;p,X,j,a]. 

This is why we may expect I e and I l to vanish simultaneously. 

What we can actually prove is the equality of the indices 1^; and \e under some ‘generic assumptions’. 

Theorem 7. Assume that {p,u) and {v,w) are periodic solutions to, respectively, (33) and (34), depend¬ 
ing smoothly on the parameters (p,X,j,a), as well as their periods S and T. With notation introduced 
earlier, ifOfj,u ^ 0, and the matrices Ce and Cl are nonsingular then we have the identity 

n(Hess0) - 2 = n(AE) - n(C£;) = n(AL) - n(CL). 

Proof. As already noticed, we have 

Ib := n(AE) - n(C£;) = n(a£;) - n(C£;) and II := n(AL) - n(CL) = n(aL) - n(CL). 
Furthermore, by Proposition 1, we have 

n(Hess6) = n(—C^;) if 0^^ > 0, n(Hess0) = n(—C^;) + 1 if 0^^ < 0, 
n(Hess0) = n(—C^) if 0 aa > 0, n(Hess0) = n(—C^) + 1 if 0 aa < 0. 
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We claim that, by Lemma 2, these relations imply 

n(Hess0) = n(-C£;) + n(a£;) - 1 = n(-CL) + n(ai) - 1, 

from which we arrive at our final formula — similarly as in the proof of Theorem 5 — by observing that 
for the 3x3, noninsingular matrices Ce and we have 

n(-CE) = 3-n(CE), n(-CL) = 3 - n(CL). 

It just remains to check that Lemma 2 does imply that 


n(aL) 

= 1 if 0AA > 0 , 

n(aL) 

= 2 if 0AA < 0 , 

(47) 

n(a£;) 

= 1 if 0^^ > 0 , 

n{aE) 

= 2 if 0^^ < 0 . 

(48) 


This a matter of adapting notation, and checking the relationship between the energy levels for the 
vector-valued profile equations and for the reduced, scalar profile equations. 

Eliminating w from the first equation in (34), we can view it as an Euler-Lagrange equation for the 
energy 

‘^{v,Vy) ■.= e{v,Vy) - - (-^-bjcr)v, 

and eliminating w from the third equation in (34), we find that the corresponding energy level is 

= A — . 

Therefore 0 aa = Ta is indeed the partial derivative of the period of u with respect to energy level, and 
since = Hess'^jw] = Hess£[u] — Lemma 2 applies and shows (47). 

We proceed similarly to prove (48). By eliminating u from the first equation in (33) we find 

^ = -A -b , 

which is the Euler-Lagrange equation associated with the energy 

'^{p,Px) :=S{p,p^) - f/{2p) - (-A-bi(T^)p, 

and by eliminating u from the third equation in (33), we see that the corresponding energy level is 

\-^[q] = p-ja. 

Therefore 0^^ = is the partial derivative of the period of p with respect to energy level, and since 
aE = Hess<f[^] — j^lp^, Lemma 2 applies and shows (48). 

□ 


5.2.2 Connexion with quasilinear KdV equations 

As already mentioned in the proof of Theorem 7, the profile equations in (34) for traveling wave solutions 
{v,w) = {v,w){y + jt) to (EKL) reduce, after elimination of w to 

E^b] -fy + p,-ja = 0, 

which is nothing but the governing ODE for the profile of traveling wave solutions v = v{x + j'^t) to 
(qKdV), and the Euler-Lagrange equation associated with 

£[v]ja - p, -f] = e{v, Vy) - + {p - jcr) v , 
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the Lagrangian defined in Section 5.1. Similarly, the velocity w can be eliminated from the abbreviated 
action integral 

= / {£{v,Vy) + — jvw +nv — (TW +X) dy 

Jo 

for (EKL). This yields 

0(M,A,j,cr) = [ {£{y,Vy) - + (y-j(j)v +^a‘^)dy = e{X- Xa'^Ja - , (49) 

Jo 

where 9 is the abbreviated action integral associated with the T-periodic traveling wave solutions v = 
y{x+j‘^t) to (qKdV), as in Section 5.1. Therefore, by the chain rule, Hess0(^, A, j, cr) can be expressed in 
terms of and Hess0 evaluated at (A — ^a^,ja — y, —j^). This will be used in §5.2.3 below to compare 
the stability criteria for (EKL) and (qKdV). 

5.2.3 Orbital stability in the Euler Korteweg system 

Let us now examine in which situation we may apply Theorem 3 to the Euler-Korteweg system. We 
consider energies as in (3), that is 

^ = P{p) + 2 "^ ip)pl ; 

or equivalently, 

£= fiv) + ■^k{v)vI , 

with 

Fip) = pf{l/p), JXX{p) = p-^K{l/p ), 

with / : (0,+oo) —)■ M and n : (0,+oo) —)■ (0,+oo) of class C^, hence also F : (0,+oo) —>■ M and 
JF : (0,+oo) —>■ (0,+oo) are of class C^. Even though many others are possible, we give now the most 
classical kind of nonlinearities that we may think of : 

• Shallow-water pressure law: F{p) = or equivalently f{v) = l/{2v), which gives p(u) = l/(2v^); 
then the EKE system is a dispersive modification of the Saint-Venant equations for shallow water 
flows, in which the ‘pressure’ term is indeed known to be of the form p{p) = ^p^; 

• NLS capillarity: JF{p) = l/(4p); then the EKE system is the fluid formulation^^ of the Non Linear 
Schrodinger equation 

(NLS) + Idl^ = MW"), 

with g{p) = F'{p). In the shallow-water case, g{p) = p, which corresponds to the cubic NLS. 

• Semilinear EKE: JF constant. 

• Semilinear EKL: k constant. 

We warn here again the reader that, as explained in §5.2.1, our notation from the abstract setting 4.1 
have to be adapted to parameters (p, A, j, a) used to define the action in (35) for both (EKL) and (EKE) 
(see in particular Table 2). 

By the connexion made in §5.2.2 between periodic waves in (EKL) and (qKdV), and the one-to-one 
correspondence between periodic waves in (EKL) and (EKE) recalled in §5.2.1, we have a good knowledge 
of situations in which (HO) is satisfied for both (EKL) and (EKE). Proposition 2 exemplifies one such 
situation. More precisely, under the assumptions of that proposition, (HO) holds true for both (EKL) 
and (EKE) for (p, A, j, a) in the preimage by 

(M,A,j,cr) 1-^ (A - - p,-f) 

the Madelung transform: ip = (px = u- 
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of the open set f2 defined in Proposition 2. Other possibilities, with various families of periodic waves, 
are considered in [BGNR14]. 

Taking (HO) for granted, let us look at our other main assumptions (H1)-(H2)-(H3). 

The abstract, nondegeneracy assumptions in (HI) become, for (EKE) 


©mm 0 ! det C^; 7^ 0 , 


and for (EKL) 

©AA 7 ^ 0 , det Ci 7 ^ 0 . 

Both of them can be reformulated in more convenient way by using Proposition 1, in which Eq. (30) 
applied to (EKE) and (EKL) gives 


©^^ detC^; = 0AA detCi = — det(Hess0). 


(50) 


(This identity may also be checked from the expressions of Cb and given in (43) and (44) respectively.) 

In particular. Proposition 1 implies that, if 0^^^ 0aa = 2^ T\ is nonzero, the matrices Cx and are 

simultaneously nonsingular, and this happens when Hess0 itself is nonsingular. In other words, (HI) 

amounts to 

©MM > det(Hess0) ^ 0 

for (EKE), and 

©AA 7 ^ 0 , det(Hess0) yf 0 

for (EKL). 



Regarding (H2), we claim that Hx := iL^(K/TZ) x L^(IR/TZ) and Hh := iL^(IR/SZ) x L^(IR/E!Z) do 
the job, respectively for (EKL) and (EKE). For the former, this comes from a straightforward addition 
to what is done in §5.1, since the Hamiltonian functional is just 


.T .T 

{v,w) / £{v,Vy)dy + I w'^dy = E[v] + ^ 

Jo Jo 


For the latter, the Hamiltonian functional is 

(/9,m)h> f ■M’{p,u,px)dx = f S’{p,pj;)dx + ^ 


\w\\h ■ 


pv? da;. 


The first part of this functional can be handled exactly as before, since £ and S’ have the same abstract 
form. By this, we mean that the mapping p i—>■ f“S(p,px)dx is on JI^(R/SZ), and for S-periodic 
P & Cl with values in (0, +oo), 

p^ {{HessS[p]p,p) + a\\p\\l 2 )^/^ 

defines an equivalent norm on iL^(IR/SZ), provided that a is large enough. Furthermore, 


^Hess.^[p,u] ^ ^ ^ ^ ^ ^ ^ = (Hess<?[p]p,p) + J {p + 2upu) dx , 

and _ _ 

^ (pu^+ 2upu)dx + a\\p\\l 2 = 

where 

A:=( ^ M \ 

“ \ H QU J 

is bounded, and positive for a large enough, its upper and lower bound being uniform when the image 
of (p,u) is in a compact set of (0, +oo) x M. This proves that (H2) holds true. 

Finally, in view of earlier work [BGDD06, BGDD07] on the Gauchy problem for (EKL) and (EKE) 
in one space dimension, we claim that (H3) is satisfied in x for s > i/2. 
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Theorem 8. Let s > 3/2. If f : (0,+oo) —)■ M and k : (0,+oo) —)■ (0,+oo) are then, for 

all T > 0, {vo,wo) € II^~^^(R/TZ) x H^(R/TZ), vg > 0, there exists T > 0 and a unique {v,w) G 
C([0, T); x iJ®(IR/TZ)) solution to (EKL) with e. = f{v) + ^K{v)vy, and {yg,wg) i—>■ iv,w) 

is continuous. 

Theorem 9. Let s > 3/2. If F : (0,+oo) — >■ M and JF : (0,+oo) — )■ (0,+oo) are then, for 

all S > 0, (po,Uo) G i7^+^(IR/SZ) x H^fR/EZ), pg > 0, there exists T > 0 and a unique {p,u) G 
C([0, T); iJ®+^(IR/SZ) X iJ*(IR/E!Z)) solution to (EKE) with S’ = F{p) + ^JF{p)p^, and {pg,ug) i-G {p,u) 
is continuous. 

We thus have all the ingredients to apply Theorem 3 to (EKL) and (EKE). 

Theorem 10. Under the assumptions of Theorem 8, if (HO) is satisfied, then for a wave profile (v, w) 
at which we have 

(Sl) Oaat^O, det(Hess0) 7 ^ 0, n(Hess0) = 2, 
the associated periodic wave solution to (EKL) is conditionally, orbitally stable in H^{R/TZ) x L^(IR/TZ). 

Proof. Our assumptions include (HO), imply (H3) by Theorem 8, and are designed to also meet (HI) 
and (H2) in the way explained earlier. Furthermore, the kernel of is spanned by (vy,w_y)^ because 
the kernel of is spanned by Vy, which follows from the first statement in Lemma 2, applied as in the 
proof of Theorem 7. Moreover, by Theorem 7 we have 

n(AL) — n(Ci) = n(Hess0) —2, 

which equals zero by assumption. This means that Theorem 3 does apply here. □ 

Theorem 11. Under the assumptions of Theorem 9, if (HO) is satisfied, then for a wave profile {p,u) 
at which we have 

(Se) det(Hess0) 7 ^ 0, n(Hess0) = 2, 

the associated periodic wave solution to (EKE) is conditionally, orbitally stable in II^{R/TZ) x L^(IR/TZ). 

Proof. It is of course very similar to that of Theorem 10. Our assumptions ensure that (H0)-(H1)-(H2)- 
(H3) are all met. The kernel of Ae is spanned by (px^tAx) because the one of is spanned by px, 
which follows from the first statement in Lemma 2, applied as in the proof of Theorem 7. The latter also 
shows that 

n(AL) — n(CL) = n(Hess0) — 2 = 0 

by assumption. So we can apply Theorem 3 here too. □ 

We may now take advantage of the connexion between (EKL) and (qKdV) pointed out in §5.2.2 to 
compare the stability criteria in Theorems 4, 8, 9. 

As shown in Appendix A, we have 

det(Hess0) = (20^6'^ - Of) (6»^^6»aa - 6*^) “ det(Hess6»), (51) 

/ 0 0 0 \ 

n(Hess0) = 1 + n(Hess6l-(26»c6»;x-6»i)/(4j26»^)J), J := 0 0 0, (52) 

Vo 0 ij 

where the derivatives of 9 are all evaluated at (A — |(t^, jct — p, —p). The fact that is a nonconstant 
periodic function of period = T implies at the same time that > 0 and 
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by the Cauchy-Schwarz inequality. Therefore, the (qKdV) stability conditions in (s) automatically imply, 
by (51), that det(Hess0) > 0 if j 7 ^ 0 and — 9\^) < 0. In this situation, we also have 

n(Hess6» - {29^9- 9\)/{‘if9^)J) = n(Hess0), 

since the involved 3x3 matrices have the same first two principal minors, and their determinants have 
the same sign. So the third condition in (s) implies n(Hess0) = 2 by (52). In fact, as soon as j 7^ 0 and 

det(Hess0) det(Hess0 — {29c9^ — 9\)/{Aj‘^9^)J) > 0, 

or equivalently — by Eq. (51) — 

det(Hess0) det(Hess0) < 0, 

Eq. (52) implies that the stability conditions in (Sl) are equivalent to those in (s), evaluated at 
(A — j(j — /i, —p)- We may summarize this in the following statements. 

Theorem 12. Under the assumptions of Theorems 5 and 10, if 


j 7 ^ 0 , 0^^ 7 ^ 0 , det(Hess0) det(Hess0) < 0 , 

then 

n(Hess0) = 1 n(Hess0) = 2 . 

If this is the case then the periodic traveling wave {v,w) = {v,w)(y + jt) is an orbitally stable solution to 
(EKL) and v = v{y + j'^t) is an orbitally stable solution to (qKdV). 


Corollary 1. Under the assumptions of Theorems 5 & 10, if j 0 and 


(s 2 ) 9^^<0, 


9\\ O^x 

9x^ Ofifi 


< 0, det(Hess0) < 0, 


holds true at (A — , ja — pL,—p), then the periodic traveling wave (v,w) = (v,w)(y + jt) is an orbitally 

stable solution to (EKL) and v = vfy + j^t) is an orbitally stable solution to (qKdV). 


6 Numerical investigation of specific examples 

In this section, we focus on the quasilinear Korteweg-de Vries equation (qKdV) and on both versions of 
the Euler-Korteweg system, (EKL) or (EKE), with an energy as in (3). In all cases, the traveling profiles 
y are governed by an equation of the form 

^K{v)y^+ W{y;X,c) = y, (53) 

where c is linked to the speed of the wave, and A, y are basically constants of integration. For (qKdV), 
notation in (53) are perfectly consistent with those used in §5.1. Indeed, c is exactly the speed of the wave, 
A is the Lagrange multiplier in the profile equation, y is the energy level resulting from the integrated 
version of this profile equation, and the potential is 

W{v;X,c) = — f(v) — + Au. 

We warn the reader once more, however, that the abstract notation for the parameters in (53) has to be 
interpreted differently for the EK system. Indeed, as should be clear from the discussion in §5.2.2, the 
profile equations for (EKL) are not readily of the form (53), which is a reduced profile equation, obtained 
only after eliminating the velocity profile w. This yields the following table of correspondence (Table 3) 
between the ‘practical’ parameters {y,X,j,a) of the EKL wave profile {y,w) introduced in §5.2.1, and 
the abstract parameters {y, A, c) in (53) — which can also be checked from the relation between the EK 
action and the qKdV action in (49). 
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abstract 


A 

c 

qKdV 


A 

c 

EKL 

A - 

jo- - M 



Table 3: Parameters in the reduced profile equation 

In fact, the computations performed in the proof of Theorem 7 also show that the reduced profile 
equation for (EKE) has an abstract form as in (53), up to adapting notation. The correspondence is 
summarized in the next table (Table 4). 


equation 

energy level 

potential W 

qKdV 


W(v) = — f{v) — -1- An 

EKL 

A-ia- 

l^{v) = - f{v) + + {ja - g)v 

EKE 

IJ--jo 

W{p) = - F{p) + ijVp -f (^CT^ - X)p 


Table 4: Energy level and potential in the reduced profile equation 
We find it convenient to keep the ‘abstract’ notation in (53) for the discussion that follows. 


6.1 Methodology 

We wish to investigate numerically the conditions for orbital, coperiodic stability that are derived in 
Section 4, and applied to (qKdV) in §5.1, and to (EKL), (EKE) in §5.2.1-5.2.2. Equation (53) is obviously 
solvable by separation of variables, which readily gives the action 


0 = 


k{v) i) dn 



K{v){fi — W{v, A, c)) dn , 


(54) 


and the period 



dv 

v/2(^- ir(v; A,c))/K(n) 


(55) 


To actually compute these integrals, we first need the points V 2 , V 3 , as denoted in Table 1, at which W 
achieves the energy level /i and in between which W is less than /x, the point V 2 being the trough of the 
wave, and V 3 , the crest. These points we can compute numerically by means of, e.g. Newton's method. 
This is what we have done, with a relative toleranee e = 10“^*^. Then we have computed numerically 


0 = 2/ du, 

J Vo 


and also 

^ ^ f'’^ _dw_ 

Jv 2 \/2{g,-W{v\\,c))/K{v) ’ 

by standard numerical integration techniques^®. Some care is needed to handle properly the square root 
singularity at endpoints. This has been done by first making the desingularizing change of variables 

V3+V2 V3- V2 . 

V = --sm uj 

2 2 


^®Here, the trapezoidal rule. We have also tried Simpson’s rule without any significant benefit. 
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(already used for instance in [BJIO, BGNR14]). Then every Newton-Cotes formulae that we tried worked 
well, and gave very close results, with a step size Aw = 10“^. We just had to preferably use an open 
Newton-Cotes formula to avoid too much sensitivity on the numerical error made on the preliminary 
computations of V 2 and V 3 . 

Finally, we have computed approximations of the second order derivatives of the action by means of 
finite differences. To avoid any confusion, the qKdV action is denoted by 9 in what follows — as in §5.1, 
and the notation 0 is reserved for the EK action. The second order derivatives of 9 and 0 are computed 
in the natural, ‘abstract’ variables, respectively and (/i, Ai, A 2 , c), up to referring to Table 2 for 

computing derivatives of the EK action with respect to the ‘concrete’ parameters (/i, A,j, cr). In both 
cases, we have used a 9-point stencil for the computation of those second order discrete derivatives, 
with a step size Aiy which turned out to be rather severely limited by the earlier steps — especially the 
numerical integration step. An additional difficulty is that the condition number of the Hessian of the 
action is most often very large. 

The tolerance e = 10“^° and the integration step size Aw = 10“"^ are basically kept constant in the 
results presented below, whereas the finite difference step Av varies from place to place. 

One important observation to validate our numerical approach is that, for some specific nonlinear¬ 
ities which correspond to completely integrable PDEs, we have access to an explicit formula for some 
characteristics of the Hessian of the action. To be more precise, we know from [BCNR13, Proposition 2] 
that the characteristic velocities of the modulated equations are those of a matrix known in terms of the 
Hessian of the action, the velocity and the period of periodic waves. For (qKdV), this matrix reads 


/ 1 0 0 

d := c 0 1 0 

\ 0 0 1 


(Hess6*) ^ 


0 0 -1 \ 
0 10 
-10 0 / 


and for (EKL), it is 


/ 1 0 0 0 
D = -7 0 1 0 0 

\ 0 0 0 1 


(Hess0) ^ 


0 

V-1 


0 0 
0 1 
1 0 
0 0 


-1 

0 

0 

0 


\ 

/ 


So, in cases when these characteristic velocities are known through ‘explicit’ formulas — which involve 
elliptic functions —, we can compare them with the eigenvalues of d and D that we obtain through our 
numerical approach. This is what we have done for the ‘standard’ Korteweg-de Vries (KdV) equation 
(p(u) = \v^, K = 1) and for the Euler-Korteweg version of the cubic NTS equation {p{v) = l/(2v^), 
J^(p) = 1/(4/9), that is, k{v) = l/(4v^)), by using the expressions of characteristic velocities collected in 
the book by Kamchatnov [KamOO] (chapter 3.5, p. 184 and chapter 5.1, p. 238) and in [Jenl4, § 2.2]. 

Remark 10. This testing against explicit formulas incidentally enabled us to find out that the numerical 
computations displayed in [BGNR14] were, to some extent, corrupted, because of inconsistent choices for 
the finite difference step Ah' and for the integration step size Aw — in fact, Av was taken too small. 


6.2 A few numerical results 

6.2.1 Benchmarks 

Let us start with the KdV case (with actually p{v) = 3v^, for convenience). On Eigure 3 hereafter, we 
have plotted the condition number of Hess0 for a reasonable range of periods T, almost going down to the 
harmonic period — that is, the small amplitude limit — as well as the relative error of its determinant 
when compared to what is expected from the explicit formulas mentioned above. As we can see, it is for 
periods close to the harmonic period that the condition number of Hess6* is the highest, and the error 
on the determinant reaches the highest values, around 0.12%. Were the step size At/ diminished, this 
instability phenomenon would worsen. 
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c 
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period 


Figure 3: Left: condition number of Hess0 as a function of the period T for KdV, with c = 60, A = —60 
kept fixed. Right: relative error on det Hess0. Finite difference step size: Ai/ = 0.005. 


Figure 4 hereafter displays, with the same parameter values as in Figure 3, the minors that encode 
the orbital stability condition in Theorem 5, namely mi = m 2 = 6 ^^6x\ — and m 3 = det Hess0. 





Figure 4: Upper left: mi = 9^^] upper right: m 2 = and lower left: m 3 = det Hess0, all of 

them as functions of the period T for KdV with c = 60, A = —60 kept fixed. Finite difference step size: 
Ai/ = 0.005. 


By looking at extreme values of those minors given in boxes, or by zooming on and eye-checking the 
parts of their curves that look closest to zero, we recover the signs that confirm orbital stability for KdV. 
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More precisely, we find that the signs of minors of Hess0 are as in the third row on Table 5 




Oxx O^x 
QX)! 


< 0 , 


det(Hess0) < 0, 


(which is (31), corresponding to the set of conditions found by Johnson in [Joh09]). 

Let us now consider the cubic NLS nonlinearity, whose numerical behavior differs. Again, on Figure 
5 below, we have plotted the condition number of Hess© and the relative error of its determinant. 




Figure 5: Left: condition number of Hess© as a function of the period for NLS, with j = 1, cr = 0, 
/X = —2.5 kept fixed. Right: relative error on det Hess©. Finite difference step size: Azx = 10“®. 


Here the most important numerical problems come from large periods, that is, when we approach the 
soliton limit. The condition number gets quite large, and the error on the determinant obviously blows 
up in this limit. Nevertheless, we can plot the minors encoding orbital stability as long as the period is 
not too large. We first plot, in Figure 6 , the four minors of Hess© encoding the stability conditions in 
(Sl) (in Theorem 10), with in particular Mi = Q\x and M 4 = det(Hess©), and then also look at the 
missing one for checking the stability conditions in (Sg) (in Theorem 11), namely ©^i^x. 
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Figure 6: The four principal minors: Mi = 0 aa (upper left); M 2 = (upper right); M 3 

(lower left); and M 4 = det(Hess0) (lower right), as a function of the period for NLS case, with j = 1, 
cr = 0, /i = —2.5 kept fixed. Finite difference step size: Az/ = 10“®. 



Figure 7: 0^^ as a function of the period with the same data as in Figure 6 . 

Again, the signs are not easily visible for periods approaching the small amplitude limit, but the 
extreme values of minors displayed in boxes show that 

0/in > 0, 0AA > 0, det(Hess0) > 0, n(Hess0) = 2, 

(the latter being a consequence of Sylvester’s rule and the fact that there are exactly two sign changes in 
the sequence of minors, since M 2 < 0 and all the others are positive), which imply that both (Sl) and 
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(Se) are satisfied. This confirms, as known from [GH07a], that cubic NLS periodic waves are orbitally 
stable, including in the mass Lagrangian coordinates of the fluid formulation of NLS. 

6.2.2 More qKdV test cases 

Let us now discuss (qKdV), with the aim of deriving new stability results in non integrable cases. We set 
again k = 1 — which means we actually concentrate on (gKdV) — and we look at a pressure law given 

by 

p{v) = e (7 + 1 )'^"'^, 7 > 2 , e = ±1 . 

The minus sign is usually referred to as the defocusing case. Its introduction is in fact irrelevant if 7 
is an even integer, because in this case the symmetry (a;, t, v) 1 —> (—x, —t, —v) drives back the minus sign 
to a plus sign in (gKdV). 

The case 7 = 2 corresponds to (KdV), while the plus sign with 7 = 8 corresponds to what is known 
as the modified KdV equation (mKdV). This is another completely integrable case, for which the results 
are qualitatively similar to those of KdV as displayed on Figures 8-9 here after. According to [AP07, 
Theorem 5.2], dnoidal waves should be orbitally stable. We chose our parameters in order to observe 
this particular family of waves (see Figure 2.1 in [Joh09]), and we find indeed the expected stability, the 
minors of Hess6* satisfying again (31). 



Figure 8: Condition number of Hess0 as a function of the period for focusing (mKdV), that is 7 = 3 and 
e = 1, with c = 1000 and A = —500 kept fixed. Finite difference step size: Av = 0.05. 
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Figure 9: Upper left: mi — 9^^] upper right: m 2 = 9^^6\x — and lower left: m 3 = det Hess0, as a 
function of the period for focusing (mKdV), with the same data as in Figure 8 . 


When going to the defocusing mKdV case (with a minus sign in the pressure and 7 = 8 ), we observe 
on Figure 11 that the intermediate minor m 2 has the opposite sign, compared to what happens for the 
focusing mKdV. However, the negative signature of Hess0 remains equal to one, and this confirms orbital 
stability for the defocusing mKdV equation, as expected from [DNll], 



Figure 10: Condition number of Hess0 as a function of the period for defocusing (mKdV), that is 7 = 3 
and e = —1, with c = —100 and A = —60 kept fixed. Finite difference step size: /S.v = 0.005. 
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Figure 11: Upper left: mi = 0^^; upper right: m2 = O^^Ow — 6\^-, and lower left: m3 = det Hess0, as a 
function of the period for defocusing (mKdV), with the same data as in Figure 10. 


Let us come back to the focusing gKdV, with this time 7 = 4. This is a non integrable case, for which 
we are not aware of any analytical result regarding the stability of periodic waves. It is only known from 
[PW92] that solitary waves are orbitally stable, which is a necessary condition for the stability of periodic 
waves of large period, by the work by Gardner [Gar97]. Numerical difficulties arise here at both ends. We 
have indeed large values of the condition number in both the small amplitude limit and the soliton limit, 
as can be seen on Figure 12. Nevertheless, there is numerical evidence for stability (Johnson’s conditions 
in (31) are satished) at least for intermediate periods, for which the numerical results are most reliable, 
see Figure 13. 



Figure 12: Gondition number of Hessd as a function of the period for (gKdV) with 7 = 4 and e = 1, and 
c = 1000, A = —500 kept fixed. Finite difference step size: Az/ = 0.005. 
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Figure 13: Upper left: mi = 0^^; upper right: m 2 = O^^iOw — 9\^-, and lower left: m 3 = det Hess0, as a 
function of the period for (gKdV) with 7 = 4 and e = 1, with the same data as in Figure 12. 


6.2.3 More EK test cases 

Let us now focus on the Euler-Korteweg system. As for the NLS case, we investigate the orbital stability 
conditions for both (EKL) and (EKE), so as to check that stability occurs at the same time in the two 
formulations. First, we consider a Boussinesq pressure law 

p(v) = V — v'^ , 

with a constant capillarity k = 1. In what follows, we take 7 = 2 , which corresponds to the good 
Boussinesq equation, as in [BJK14, § 4.2.1]. 

Numerical results are displayed on Figures 14-15. We observe a transition at period Tq ~ 3.68, where 
det Hess© changes sign (and n(Hess0) passes from 2 to 3). For periods smaller than Tg, we see that 
we are in the range of application of Theorems 10 & 11, which imply orbital stability for both (EKL) 
and (EKE). For periods larger than Tg, since det Hess© < 0, Theorem 1 implies spectral instability. We 
checked that the zone around Tg where the condition number becomes very high does not depend on 
our choice of discretization steps — unlike what happens in the KdV/NLS cases for periods approaching 
the harmonic one or going to infinity. The transition at Tg is not a numerical artifact. Our conclusions 
are thus consistent with those in [HSS14], where it is pointed out that, at a fixed velocity, there exists a 
maximal period for the wave to be stable. 
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Figure 14: Left: condition number of Hess© as a function of the period for EKL with Boussinesq pressure 
law with j = 2, j = —0.1, cr = 0, /r = —2 kept fixed. Right: Finite difference step size: 

Ai/ = 0.5 • lO"-*. 






Figure 15: The four principal minors: Mi = Q\x (upper left); M 2 = ©^^©aa — ©a^x (upper right); M 3 
(lower left); and M 4 = det(Hess©) (lower right), for EKL with Boussinesq pressure law with 7 = 2 and 
the same data as on Figure 14. 


Let us now consider a last case. The perfect gas pressure law 

p{v) = l/ 2 w , 

with a constant capillarity J^{p) = 1 , or equivalently, k(v) = 1/v^. 
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Figure 16: Left: condition number of Hess© as a function of the period for EK with the perfect gas 
pressure law, with j = —1, u = 0, /r = —2.5 kept fixed. Right: Finite difference step size: 

Ai/ = 0.5 • lO"-*. 






Figure 17: The four principal minors: Mi = ©aa (upper left); M 2 = ©^^©aa — ©a^x (upper right); M 3 
(lower left); and M 4 = det(Hess©) (lower right), for EKL with the perfect gas pressure law and the same 
data as in Figure 16. 


We observe that (Sl) and (Se) are satisfied, so that both Theorems 10 & 11 apply. The tested 
periodic waves in the Euler-Korteweg system with the perfect gas pressure law are orbitally stable in 
both mass Lagrangian coordinates and Eulerian coordinates. 


54 




















Appendix 


A Algebraic computations regarding 0 and 6 

Recall from (49) that 

= ^'(A- 5 cr^,jcr-/x,-j^), 

and that the subscripts fj,, A, c will denote, when attached to 0, partial derivatives of 9 with respect to 
its first, second, and third variable respectively. Therefore, by applying twice the chain rule we obtain 

f ^ ^\fi ^jfJ' ^afi \ 

HpqqP) = 

©AO ©A, ©-. 

\ ©//(7 ©A(7 ©jct ©crcT j 


/ 9\x 

—9xyL 

-ct6»aa + 2j6»cA 

-j9x\ + <j9xyi 



-2j9c^ + (J9xf, 

-h j0xt^ 



'^f9cc + <j‘^9xx 
-4:ja9cx - ‘2-9c 

jcr9xx - 2j‘^9cx + 2ja9cfj. 

-<j'^9xyi + 9x 

V 



f^0xx + 

-2ja9xfi - 9f, 


where the entries of Hess© are evaluated at (^, A, j, cr), while the derivatives of 9 are meant to be evaluated 
at (A— , j(T — fi, —j^) — all along this section^®—, and we have omitted to write the symmetric entries 

under the diagonal in order to save some space. 

By elementary manipulations on columns and rows — as in Proposition 1 — we thus find a matrix 
Q € SL 4 (IR) such that 


Q(Hess0)QT 


9xx 

-9xn 

2j9,x 

0 ' 

—9xn 


-2j9cf^ 

0 

2j9cx 

-2j9cf^ 

4f9,, - 29, 

9x 

0 

0 

9x 

-0, , 


Assuming that j 0 — and recalling that T = > 0 — we can perform further manipulations, and 

thus find R e GL 4 (IR) such that detR = l/(2j) and 


R(Hess0)R'r 


/ 9xx 

9xxl 

9cx 


9xfj. 

e,. 

9cyL 

0 

9cx 

dcfl 

9cc - V 

0 

V 0 

0 

0 

/ 


,7 := {29,9^ - 91)/{4jX) . 


(56) 


Observe that the upper 3x3 block in the right-hand side is Hess0 — rjj, where 

/ 0 0 0 \ 

J := 0 0 0 . 

Vo 0 ij 


In particular, (56) yields 

-4t det(Hess0) = —det(Hess0 — rjj ), 

which equivalently reads 

det(Hess0) = (29^9^ - 9l) {9^f,9x\ - 9l^) - Af 9^ det(Hess6>). 

^®In particular, the reader should keep in mind that with our conventions, 

= ®AA • 
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Since 0^ > 0, (56) also gives 

n(Hess0) = 1 + n(Hess6l — ryj). 

This proves Eqs. (51) and (52). 

We can also draw some sign tables that yield the signatures of Hessd, C^, and C^; by inspection of 
their minors. 

Let us first consider C^, written as in (44), 


Cr = -- 



/ 

^Xfi 


^ fl(T 



\ 



Qtj.X 0AA 


0/xA Baa 


0;iA 0AA 


1 


^afj, 0A/^ 


Bcrcr B^(7 


0(Tj 0Aj 


0AA 


0o'A 0AA 


Bcta Baa 


0<tA 0AA 




0j> 0A;i 


Bjcr Ba(7 





1 

0jA 0AA 


Bja Baa 


0jA 0AA 

/ 


which equivalently reads 


( 

9xx —9xti 


0 0 


2j9xc —2j9fic 

\ 


—9fix 9^^ 


^/j.A ^MM 


—9^x 9^i^ 



Cr =- 


0 


AA 


0 -9x^ 


0 

1 


9x 

-2j9^c 

0 9^^ 


0 9^^ 


0 




2 j6*cA 


9x 0 


Af9,,-29, -2j9^, 


\ 





^MM 

/ 


1 

0AA 


2 j 


^AA ^Ap 
0 

ScX O^x 


2 j 


^Ac ^[ic 
Sxfj. d/ifi 




9x9u 




9x9f^fj. 

9 cc 9 

9cfi 9f_ipi 


- 29,9^^ 


This is to be compared with the constraint matrix for (qKdV), which is given by 


f 


V 


9xx 9f^x 
9xfi 9^^ 

9cx 9^x 

9cfi 9fj,^ 


9xc 9 
9xfj. 9^^ 


9cfj, 9 


flC 

MM 




J 


Recalling that 0 aa = 9f^f^, we see that the first upper left minors of and c are both equal to 


Ai :=- 

The second upper left minor of is 

A 2 := —9^ 


1 

9„n 


9xx 9fj,x 

9xfi 9^fi 


=: 6. . 


9xx 9fj,x 
9xiJ, 9 


Hfj. 


= 9a 9an ^1 , 


82 ■= det c = 


; minor of c 

is 





9xx 9f^x 


^Xc ^ 



1 


9x^1 9^^ 


^A/j- 


1 



9cx 9^x 


^cc ^ flC 


9^,^, 



9cfi 9fj,fj, 
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As to the third and last hrst upper left minor of C^, it reads 


A 3 det Cl = 




that is, 


Qxx 

Opix 

Oxfi 

dpipi 

OcX 

Opix 

Ocfl 

Opifi 

A 3 

= 4fo 


^Xc ^ 

Sxfj. 0^^ 

^ cc G pic 
^cpi ^ pipi 

62 + {29(.0n — 



Ol)Si. 


Oxx 

(^piX 

Sxpi 

Opipi 


Since 0^ = T > 0 and 29cdfi — 9\>Q (by Cauchy-Schwarz), we can compare n(Ci) and n(c) by drawing 
the following table, which relies on Sylvester’s invariance theorem^^, and where +/— signs stand for 
positive/negative values: 


dfipi 

O^fiOxx 

det(Hess 0 ) 

n(Hess 0 ) 

< 5 i 

^2 

n(c) 

Ai 

A2 

A3 

n(CL) 

det(Hess 0 ) 

+ 

+ 

+ 

0 

- 

+ 

2 

- 

- 

? 

> 1 

? 

+ 

+ 

- 

1 

- 

- 

1 

- 

- 

- 

1 

+ 

+ 

- 

- 

1 

+ 

- 

1 

+ 

+ 

? 

> 0 

? 

+ 

- 

+ 

2 

+ 

+ 

0 

+ 

+ 

+ 

0 

- 

- 

+ 

+ 

2 

+ 

- 

1 

+ 

- 

? 

> 1 

? 

- 

+ 

- 

3 

+ 

+ 

0 

+ 

- 

+ 

2 

+ 

- 

- 

- 

1 

- 

+ 

2 

- 

+ 

? 

> 2 

? 

- 

- 

+ 

2 

- 

- 

1 

- 

+ 

- 

3 

- 


Table 5: Possible values of signs of minors and negative signatures regarding Cl 


We have used here above the additional observation (already made in [BGNR13, Remark 2]) that 



^AA ^/iA 


^Ac ^/i.c 


^A/j, 


^A/i. 


^cA ^a^A 


^CC ^ (1C 


^c/i 


^C(l ^ 


Opipi det(Hess0). 


Similarly as for Cl, we can express 



/ 

Baa B^a 




Ba(7 ^ (la 

\ 



Ba// B^^ 


0A^ ^pifi 


Ba/j, B^^ 


1 


Bja B^a 




Bj{j Byj,CT 


Qfipi 


Bj^ B^^ 




Bj^ B^^ 




Bo-A Byj,A 


0(Ti 0/ri 


B(7cr B^ct 



V 

B(7/i B^p 


0cr/r ^ fifi 


Bff/j, B^^ 

/ 


^^This theorem implies that the negative signature of an n X n symmetric matrix is the number of sign changes in the 
sequence (1, Ai,..., A„), where Aj. denotes the matrix’ fc-th upper left minor. 
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in terms the derivatives of 9 as 


1 

^/iA 


1 

to 

2j0xc 


—^A^ ^AA 


— &Xf, 

Oxx 


Cr = -- 


-2j6*^c -0Xe 


4f9„ - 20, 

2j9xc 

“^jdxc 0XX 


2j0,x 

0XX 


9\\ 


We thus find that the upper left minors of C^; are 

1 


— 

9x 

2j0ac 


0 

Oxx 

0 

9xx 


Asp := — 


AA 


9xx 9fj,x 
Ox^ " 




9aa 


A 


E,2 ■ = 


^E ,3 = det Ce = 


1 

9xx 

(Oxx) 

1 




^A/j, 


0 Oxc 


4j2 


^fiX ^XX 


Ofix Oxx 


{OxxY 


^cfi ^Xfi 


Occ Oxc 




OcX Oxx 


OcX Oxx 



20 , 

9xx 


011 . 1 ,. 0 . 


9^a 


Xfi 

0XX 


(4j^ det Hess0 + 20,0^^, ), 


Ofifi Oxfi 


Ofxc Oxc 

Ofix Oxx 


O^x Oxx 

Ocfi Oxfi 


Occ Oxc 

OcX Oxx 


OcX Oxx 


- — (20,9^ - 0l) 

PAA 


= V— (4j det Hess6» + {20,0^ - 0i)9^^, 6i ). 


9 


XX 



o 

1 


Ox 2j0^, 


O 

1 


\ 

0 Oxx 


0 Oxx 


0 0XX 

) 




A/2 

AA 


Furthermore, since 0, = /g > 0, we can complement Table 5 with a similar table in terms of the 
negative signature of Ce- 


Oxx 

OxxO^^ O'x^ 

det(Hessd) 

n(Hess 0 ) 

A_e,i 

l^E ,2 

A_e ,3 

^{Ce) 

det(Hess 0 ) 

+ 

+ 

+ 

0 

- 

? 

7 

> 1 

7 

+ 

+ 

- 

1 

- 

- 

- 

1 

+ 

+ 

- 

- 

1 

+ 

7 

7 

> 0 

7 

+ 

- 

+ 

2 

+ 

+ 

+ 

0 

- 

- 

+ 

+ 

2 

+ 

7 

7 

> 0 

7 

- 

+ 

- 

3 

+ 

+ 

+ 

0 

+ 

- 

- 

- 

1 

- 

7 

7 

> 2 

7 

- 

- 

+ 

2 

- 

- 

- 

1 

- 


Table 6: Possible values of signs of minors and negative signatures regarding Ce 

B Proof of Lemma 2 

For convenience of the reader, let us recall the statement of this lemma. 

Lemma. Assume that k : / —)■ (0, +oo) and W \ I ^ M. are on some open interval I and such that 
the Euler-Lagrange equation E‘W[v] = 0 associated with the energy 

^ : {v,Vx) \n{v)v1 + W{v) 
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admits a family of periodic solutions v taking values in I, parametrized by the energy level p, = for 

p G J, another open interval. If we denote by T the period of v, and assume that its derivative with 
respect to p, does not vanish, then the self-adjoint differential operator a := Hess'^[u] has the following 
properties: 

• the kernel of a on L^(IR/TZ) is the line spanned by v„.; 

• the negative signature n(a) of a is given by the following rule: 

* if'^fj,>0 then n(a) = 1, 

* < 0 then n(a) = 2. 

This lemma belongs to the classical periodic Floquet/Sturm-Liouville theory that may be found in 
[MW79, Part I-Chapter II], [RS78, Chapter Xlll-Section 16], [Tesl2, Part I-Section 5.6] or partially in 
[Nev09]. Some complements in [Nev09] are also used there and in [NN14] to address the question of 
co-periodic stability. 

Proof. The first point follows from an elementary, ODE argument. Indeed, the equation an = 0 is a 
second order ODE, of which we already know two independent solutions, namely v,,. and v^, obtained by 
differentiating the Euler-Lagrange equation E'^[y] = 0. In addition, by differentiating the null function 
v(- + T) — v with respect to p, we see that + T) — Since is not identically zero, and 

7 ^ 0 by assumption, cannot be T-periodic, nor can be any linear combination of and v,,., since 
the latter is T-periodic. Therefore, the set of T-periodic solutions of an = 0 is the line spanned by n,,,. 

The second point involves some classical spectral theory of self-adjoint operators in general, and of 
Sturm-Liouville operators in particular, as well as a technical argument based on a computation that is 
hardly ever done in details. The general theory says — see [Tesl2, Theorem 5.37] — that, since is in 
the kernel of a and vanishes only twice on one period, 0 is either the second or the third eigenvalue of 
the Sturm-Liouville operator 

a = -d„:Kdx + q, K := k{v) , q := P {v)v_„f). 

So we are left with showing that T^ > 0 corresponds to the first case, and T^ < 0 to the second. This 
will follow from a computation based on the discriminant associated with the eigenvalue equation 

av = rv. 

This discriminant is the function T = T(r) defined as 

T(r)=tr (f(T;r)), 

where f(-; r) denotes the fundamental solution of the eigenvalue equation av = rv viewed as the first-order 
system 

(” = 7 /*'', (57) 

w = [q — r) V . 

(We warn the reader that f and w are local notations, which do not have the same meaning as in the 
main part of the paper.) The reason why T plays a prominent role here is that the matrix value function 
f has a determinant constantly equal to one — this follows from Liouville’s formula, since the matrix in 
System (57) is traceless — so that f at x = T admits 1 as an eigenvalue if and only if T(r) = 2. In 
addition, this is exactly the condition under which (57) has a T-periodic solution. Therefore, T(r) = 2 is 
an iff condition for r to be an eigenvalue of a on L^(IR/TZ). In particular, we already know that 0 is an 
eigenvalue of a, so T(0) = 2. 

The connection of T with the ordering of eigenvalues arises through its derivative. Indeed, the sign 
of T'(0) tells us if 0 is the second or the third eigenvalue of a on L^(IR/TZ). This is a consequence of a 
deeper knowledge on the spectrum of a on L^(IR). By Floquet-Bloch decomposition — see point (d) in 
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[RS78, Theorem XIII-85] combined with [RS78, Theorems XIII-88 & XIII-89] for the semi-linear case, or 
[RodlS, p.30-31] for a general argument —, 

o'L2(K)(a) = [J , 

i/G[0,27r) 

where a.^, is defined by the same formula as a, that is := —dxKdx + q, but on the domain 

{v G H'^{R/TZ) ; u(T) = e*"^ t>(0)} . 

By a similar reasoning as for a, we can see that the eigenvalues of a^ are those r such that e®'^ is an 
eigenvalue of f(T;r) and this is equivalent to T(r) = 2 cosj2 . We thus see that CT 2 , 2 (R)(a) is the union of 
closed intervals in which T achieves values in [—2,2] (see Figure 18 below). 



Figure 18: Oscillations of discriminant T(r) (blue curve), spectrum of ap (abscissas of green bullets), and 
(part of) spectrum of a (red intervals). 


In particular, the eigenvalues of a on L^(M/TZ) — that is, of the operator ap — are found at the 
intersection of this graph with the line T = 2. In addition, it is known — see [Tesl2, Theorem 5.33] — 
that the least eigenvalue of ap, say ri, is in fact the lower bound for the spectrum of all the operators a^. 
By the mean value theorem, this implies T(r) > 2 for all r < ri. Paying closer attention to the variations 
of T, we thus infer that 0 is the second eigenvalue of ap if T'(0) > 0, and the third one if T'(0) < 0. 

It remains to make the connection between the signs of r'(0) and T^. This will follow from a fairly 
general computation, involving any two independent solutions (ui,r(;i)^ and (^ 2 ,^ 2 ) of (57), smoothly 
parametrized by r and such that 


Vi V2 
Wi W2 


(Again, this is possible since the matrix in System (57) is traceless.) The solution operator s{x,s;r), 
related to the fundamental solution by f (x; r) = s{x, 0; r), is thus given by 


v(a:)Aw(s) -v(x) A v(s) \ 

^ y w(a;) A w(s) —w{x) A v(s) J ’ 

where we have introduced, for convenience, the vector-valued functions 
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and used the notation A as a shortcut for the determinant in Note in particular that the requirement 
made on the determinant of (vi,wi)^ and (u 2 ,W 2 )^ is equivalent to 


V A w = 1. 


Furthermore, we have by definition 

T(r) = tr (s(T, 0; r)) = v(T; r) A w(0; r) + v(0; r) A w(T; r). 

Without loss of generality, we may assume that v(0; r) = v(0; 0) and w(0; r) = w(0; 0) for all r, so that 

T'(0) = Vj.(T; r) A w(0; 0) + v(0; 0) A Wr(T; r), 

where the subscript r stands for a partial derivative with respect to r. We claim that 

sign(T'(0)) = sign(w(T; 0) A w(0; 0)). (58) 

The proof of (58) proposed here may be thought of as a slight variation on the one relying on For¬ 
mula (2.17) in [MW79]. 

In order to prove (58), let us first observe that, for any solution (v, wY to (57), by differentiation with 
respect to r, {vr,WrY solves the system with source term 

j Vr = Wr/K , 

( Wr = {q — r) Vr — V , 

and is therefore given by the Duhamel formula associated with the solution operator s(a;, s; r). This gives 

( Vr{x) = fj(v(x) A v(s))u(s) ds 
[ Wr{x) = fg (w(x) A v(s))v(s) ds 

where we have omitted to write the dependence on r for simplicity, hence 

T^(0) = f (v(T) A v(s))(v(s) A w(0))-I-(v(0) A v(s))(w(T) A v(s)) ds , 

Jo 

where all functions are evaluated at r = 0. The idea is now to find some cancellations in the integrand 
here above, by using the known constraint v A w = 1, together with the fact that T(0) = 2, that is, 

v(T) A w(0) + v(0) A w(T) = 2, (59) 

where again all functions are evaluated at r = 0. By elementary algebra, these relations imply that 

(v(T) — v(0)) A (w(T) — w(0)) = 0 . 

At least one of these two vectors v(T)—v(0) and w(T) — w(0) is nonzero, because the solutions (vi,wi)^ 
and {v 2 ,W 2 Y of (57) cannot be both T-periodic (recall that the kernel of a is one-dimensional). 

Let us assume that u := v(T) — v(0) is nonzero. Then it is necessarily independent from v(0) — 
otherwise, there should exist a real number a such that v(T) = av(0), and (59) would imply a-f 1/a = 2, 
hence a = 1 and u = 0. Then we may decompose 

w(T) = w(0) -I- wu, v(s) = a(s)v(0) -I- /3(s)u, Vs . 

We substitute these expressions in T'(0), and by using the relation 

u A w(0) -I- to v(0) A u = 0 , 
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which comes from v(T) A w(T) = 1 and v(0) A w(0) = 1, we find that all but one term cancel out, in 
such a way that 

T^(0) = f a(s)^ u A v(0) ds = v(T) A v(0) f a(s)^ds. 

Jo Jo 

This also gives 

w^T'(O) = w(T) A w(0) f a(s)^ds. 

Jo 

Observe in addition that a is not identically zero if we set a(0) = 1 (and /3(0) = 0), so that a(s)^ ds > 

0 . 

The case v(T) = v(0) is even easier to deal with, because then 

T^(0) = f (v(0) A v(s))(w(T) — w(0)) A v(s) ds . 

Jo 

If we decompose v(s) = a(s)v(0) + /3(s)w(0) with q;(0) = 1 and /3(0) = 0, which is possible because 
v(0)Aw(0) = 1 so that v(0) and w(0) are not colinear, we can also write, using that v(0)A(w(T)— w(0)) = 
0 — which comes from the difference between v(0) A w(T) = v(T) A w(T) = 1 and v(0) A w(0) = 1, 

'T 

T'(0) = w(T) A w(0) f /?(s)^ds . 

Jo 

We observe that /3(s) = v(0) A v(s) is not identically zero, otherwise we would have 

0 = v(0) A v'(s) = v(0) A w(s)/itr(s) 

by the first equation in (57) applied to both components of v and this would contradict v(0) A w(0) = 1. 

We thus have (58) as soon as w(T) A w(0) ^ 0, and if w(T) A w(0) = 0 we have sign(T'(0)) = 
sign(v(T) A v(0)). We are in fact going to apply (58) in the hrst case. 

Let us now apply (58) to 

ui(-;0)=u,,, wi{-\0)=Kv^^, v2{-]Q) = v^, W2{-;0) = K . 

By differentiating the equation L'^ [v] = fi with respect to fi, we see that 


- K 


= 1 , 


which exactly means that v A w = 1, and 


w(T;0)Aw(0;0) = K(0)2^,J0)(u^,,(0)-u^_,(T)) = 

by using ~ Lp,a:(0) = fj, V.xxi^)^ which comes from the differentiation with respect to /r of 

tfa:~ La;(0) = 0. Since K does not vanish and VxxW ^ ^ — provided that v is chosen such that 
v,j,(0) = 0, this eventually shows that 


sign(T'(0)) = sign(T^). 
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